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ABSTRACT 


i  i  i 


Based  on  the  time-domain  first  order  correction  to  the  physical 
optics  current  approximation,  a  relationship  between  the  phase  factors 
of  the  polarimetric  scattering  matrix  elements  and  the  principal 
curvatures  at  the  specular  point  of  a  scatterer  is  established. 

The  above  phase-curvature  relationship  is  tested  by  applying 
it  to  theoretical  as  well  as  experimental  backscatter ing  data 
obtained  for  a  prolate  spheroidal  scatterer.  The  results  of  these 
tests  not  only  determine  the  acceptabi 1 i ty  of  the  phase-curvature 
relationship,  they  also  point  out  the  range  of  kb  values  over  which 
the  first  order  correction  to  the  physical  optics  currents  is  valid, 
and  which  serves  as  a  compromise  range  between  the  high  frequency 
condition  required  by  the  curvature  recovery  model  and  the  drawback 
to  lower  frequencies  required  to  prevent  critical  magnification  of 
measurement  errors. 

Another  curvature  recovery  equation  is  derived  by  transforming 
the  linear  polarization  basis  to  the  circular  polarization  basis. 

The  curvature  recovery  model  is  proven  to  satisfy  the  image  reconstruction 
identities  of  invariant  transformat  ion.  A  scattering  ratio  is  defined 
and  its  behavior  is  investigated  at  high  frequencies.  Its  plots  on  the 
complex  plane  provides  a  simple  way  to  help  judge  the  accuracy  of 
polarimetric  scattering  measurement,  regardless  of  whether  a  linear  or 
a  circular  polarization  basis  is  used. 
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CHAPTER  I 

INTRODUCTION  AND  LITERATURE  REVIEW 

1 . 1  I  nt  roduct i on 

When  a  conducting  object  is  illuminated  by  electromagnetic 
radiation,  in  general,  radiation  is  scattered  in  all  directions 
by  the  object.  The  problem  of  direct  scattering  is  that  of  determining 
the  scattered  field  in  all  directions  when  the  properties  of  the 
incident  field  as  well  as  those  of  the  object  are  known.  The  problem 
of  far-field  inverse  scattering  is  defined  as  that  of  extracting 
the  geometrical  properties  and/or  reconstruct! ng  the  shape  of  the 
scatterer  under  interrogation,  given  the  incident  field  within  the 
neighborhood  of  the  scatterer  and  the  scattered  far  field.  This  problem 
is  fundamental  to  problems  such  as  radar  target  classification, 
discrimination  and  identification  in  remote  sensing  and  surveillance 
(Boerner,  1978,  1980). 

It  has  been  demonstrated  that  an  electromagnetic  scatterer 
acts  as  a  sensitive  polarization  transformer,  depending  on  its  profile. 
Thus,  depolarization  effects  must  be  taken  into  account  in  this 
problem  of  vector  nature.  Sinclair  ( 1 948 )  introduced  the  scattering  matrix 
of  a  radar  target  for  complete  depolarization  character i zat ion.  The 
question  then  arises  as  to  what  geometrical  properties  of  the  surface 
profile  of  the  target  may  be  extracted  from  the  scattering  matrix 
of  the  target.  A  study  of  differential  geometry  reveals  that  for  a 
smooth,  convex  shape,  there  exists  a  pair  of  lines  of  principal 
curvatures  orthogonal  to  and  intersecting  with  each  other  at  any 
specular  point  on  the  scatterer's  surface.  Since  the  principal 
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curvatures  at  a  point  totally  determine  any  other  normal  curvatures 
at  the  same  point,  these  two  principal  curvatures  serve  as  a  complete 
curvature  characterization  at  the  specular  point  of  the  scatterer. 

It  will  be  shown  in  this  thesis  that  if  the  polarization  of  the 
incident  magnetic  field  is  in  one  of  the  directions  of  the  principal 
curvatures  at  the  specular  point  of  interest,  then  the  cross-polarized 
backscattered  returns  vanish.  This  kind  of  null  polarizations  generates 
the  idea  that  the  principal  curvatures  are  possibly  related  to  the 
non-zero  co-polarized  backscattered  returns,  and  it  also  generated  many 
recent  studies  in  radar  polarimetry.  This  specific  degeneracy  of  target 
depolarizat ion  phenomenology  provides  the  initial  motivation  for  writing 
this  thesis. 

1.2  Literature  Review 

Inverse  techniques  cover  a  vast  amount  of  literature  and  they 
have  been  developed  in  many  otherwise  diverse  fields  of  physical 
sciences  where  the  characteristics  of  a  medium  are  estimated  from 
experimental  data,  obtained  from  measurements  made  usually  at  a 
distance  from  the  medium,  utilizing  the  laws  th^t  relate  these 
characteristics  to  the  experimental  data  in  a  given  situation.  The 
inverse  problem  of  electromagnetic  scattering  has  not  been  solved 
in  the  general  case.  Stringent  requirements  are  often  needed  to  be 
kept  on  the  shapes  of  the  scatterer  to  be  recovered  and  on  the 
operating  frequency  range.  Moreover,  the  existing  solutions  generally 
demand  an  exhaustive  amount  of  input  information  of  as  many  aspects 
and  frequencies  as  possible.  There  are  many  approaches  for  obtaining 
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approximate  solution  to  the  electromagnetic  inverse  scattering  problem. 

The  following  only  highlights  the  ramp  response  approach  and  other 
specific  approaches  which  have  been  proven  successful  and  potentially 
promising  for  future  research,  and  which  are  directly  relevant  to  this 
thesis.  An  excellent  and  quite  complete  review  of  electromagnetic  inverse 
problems  is  given  in  Boerner  [21]. 

Kennaugh,  Cosgriff  and  Moffatt  have  introduced  the  use  of  the 
impulse  response  in  electromagnetic  scattering  problems  [10,12].  Using 
the  physical  optics  approximation,  it  was  shown  that  the  impulse  response 
of  a  smooth  conducting  object  is  directly  proportional  to  the  second 
derivative  of  the  cross-sectional  area  of  the  scatterer.  This  remarkable 
high  frequency  inverse  scattering  identity  is  known  as  the  Kennaugh-Cosgr i f f 
formula  [10].  The  area  profile  can  thus  be  recovered  from  the  ramp  response, 
since  the  ramp  response  is  then  directly  proportional  to  the  cross-sectional 
area  rather  than  its  derivative  (Kennaugh  and  Moffatt,  1965).  Young  [22] 
used  the  ramp  resDonse  synthesized  from  complex  scattering  data  at  ten 
harmonically  related  frequencies  in  the  target’s  low  resonance  regime  to 
estimate  the  area  function,  from  which  a  "likely  image"  of  the  target  was 
generated  at  three  orthogonal  look  angles  using  his  approximate  limiting 
surface  technique.  The  images  obtained  are  decent,  but  in  general  they  are 
not  uniquely  specified  since  more  than  one  shape  satisfies  any  three 
look-angle  area  function  set. 

A  more  systematic  approach  was  carried  out  by  Das  and  Boerner,  who 
showed  that  the  reconstruction  of  a  smooth  conducting  target,  convex  in 
shape,  can  be  considered  as  a  two-step  process:  (i)  an  electromagnetic 
step  of  obtaining  suitable  radar  measurables  from  which  one  can  make 
an  estimation  of  a  geometrical  function  of  the  conducting  scatterer;  and 
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(ii)a  geometrical  step  of  reconstructing  the  shape  and  size  of  the 
object  from  the  knowledge  of  the  geometrical  function  estimated  in 
the  first  step.  With  the  Radon  transform  concept,  they  showed  that 
the  inverse  three-dimensional  Radon  transform  of  the  cross-sectional 
area  of  a  scatterer  is  simply  the  target's  characteristic  function, 
which  is  a  complete  spec i f icat ion  of  shape.  The  area  function  can, 
as  in  the  previous  approach,  be  estimated  by  the  ramp  response  method 
(Boerner,  1980).  Thus,  the  electromagnetic  inverse  scattering  problem  can 
be  formulated  as  the  classical  Radon  problem.  Moreover,  they  also 
indicated  that  the  classical  Radon  problem  is  intimately  related  to 
the  problem  of  reconstruction  from  projections  which  has  long  been 
investigated  and  applied  in  diverse  fields,  particularly  in  Computer- 
Assisted  Tomography.  Thus,  the  many  reconstruction  techniques  and 
algorithms  well  developed  in  other  fields  can  be  applied  in  radar 
target  imaging  as  well. 

A  solution  for  the  inverse  scattering  problem  using  the  space- 
time  Integral  equation  was  reported  by  Bennett  et  al  [23].  The 
technique  developed  was  iterative  and  restricted  to  the  class  of 
rotational ly  symmetric  conducting  targets.  In  this  approach,  the 
Inverse  problem  is  formulated  as  an  inversion  of  the  space-time 
integral  equation.  The  shapes  reconstructed  are  excellent  for  this 
simple  class  of  shapes,  but  no  extension  to  more  general  shapes  has 
been  made  recently. 

An  important  aspect  of  electromagnetic  inverse  scattering  is 
to  incorporate  the  problem  with  utilization  of  polarization.  A 


5 


monostatic  inverse  scattering  model  based  on  polarization  utilization 
was  developed  by  Chaudhuri  and  Boerner  [2,17.2**].  In  brief,  approximation 
was  made  of  high  frequency  scattered  fields  and  an  equivalent  ellipsoidal 
model  was  developed.  The  utilization  of  the  space-time  integral  equation 
and  Minkowski's  theory  lead  to  a  system  of  equations  for  the  recovery 
of  the  surface  profile.  The  numerical  technique  is  iterative,  and  when 
the  curvature  difference  at  the  specular  point  approaches  to  zero,  the 
recovery  of  that  particular  point  is  not  possible;  the  system  of  the 
recovery  equations  also  becomes  ill-conditioned  even  if  the  curvature 
d i f ference  i s  sma I  I . 

Based  on  the  first  order  polarization  correction  to  physical 
optics  [8],  Ho  graphically  reconstructed  the  shape  of  a  sphere-capped 
cylinder  with  polarization  correction  incorporated  into  the  Radon 
transform  approach  [25,26],  The  results  show  that  the  quality  of 
images  are  significantly  improved  with  polarization  correction.  It 
must  be  noted  here  that  Ho  took  advantage  of  the  plane  symmetry  of 
the  sphere-capped  cylinder,  and  reduced  the  Radon  transform  by  one 
dimension.  In  this  class  of  objects,  the  two-dimensional  inverse 
transform  of  the  area  function  normal  to  the  equatorial  plane  gives 
the  width  perpendicular  to  the  plane.  The  width  over  this  plane  of 
symmetry  is  actually  a  complete  specification  of  the  shape.  Since 
in  the  two-dimensional  case,  the  Radon  transform  becomes  the 
projection,  thus  the  filtered  back  projection  algorithm  was  directly 
borrowed  from  the  theory  of  reconstruction  from  projections  (Shepp 
and  Logan,  107*0  . 

lr,  this  thesis,  a  high  frequency  inverse  scattering  model  is 
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developed  for  the  recovery  of  the  specular  point  curvature  from 
polarimetric  scattering  data.  Not  only  does  the  model  show  that  the 
specular  geometry  can  be  directly  extracted  from  polarimetric 
data,  but  it  also  contributes  in  viewing  the  electromagnetic 
inverse  scattering  problem  as  one  of  reconstruction  from 
curvatures  in  differential  geometry.  In  the  introduction  of  the 
next  chapter,  the  main  objective  of  this  thesis  will  be  clearly 
spec! f ied. 
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CHAPTER  I  I 

THE  SPACE-TIME  INTEGRAL  EQUATION  AND  FIRSl  ORDER  CORRECTION 
TO  PHYSICAL  OPTICS 

2. 1  I nt  roduct i on 

In  recent  years,  due  to  advances  made  in  radar  technology,  it 
has  become  possible  to  measure  the  complete  relative  phase  scattering 
matrix  of  an  object  reliably  [1,2,3].  Thus  the  utilization  of  these 
polarimetric  scattering  data  in  radar  target  identification  or 
discrimination  and  in  other  inverse  scattering  applications  has 
become  of  considerable  interest  in  current  theoretical  and  experimental 
research  efforts  [ 1  - 7 ] •  The  main  objective  of  this  thesis  is  to 
investigate  the  information  content  of  the  scattering  matrix  [S],  on 
the  shape,  size,  curvature,  etc.  of  a  scatterer,  when  it  is  given 
for  the  monostatic  case  and  in  the  high  frequency  region  (i.e.  the 
wavelengch  of  the  interrogating  signal  is  small  compared  to  the 
object  characteristic  dimension). 

The  scattering  matrix  [S],  which  manifests  total  polarization 
information  for  a  fixed  frequency  and  a  given  aspect,  is  comprised 
of  four  measurable  complex  elements  (four  magnitudes  and  four  phases). 
It  will  be  shown  here  that  the  difference  in  suitable  phase  terms 
in  this  matrix,  under  the  high  frequency  i nterrogat i on  conditions, 
can  lead  to  the  recovery  of  the  difference  in  principal  curvatures 
at  the  specular  point,  from  a  given  general  [S1  matrix.  This  procedure 
avoids  unitary  transformat  ions  used  in  the  pursuit  of  cross-polariza¬ 
tion  nulls  required  in  certain  radar  target  ident i f i cat  ion  techniques 
[ 3~ 5 J •  The  underlying  concept  used  to  achieve  the  above  results  is 
based  on  Bennett's  first  order  far-field  polarization  correction  [8,9] 


to  the  Kennaugh-Cosgr i f f ' s  physical  optics  formula  for  the  electro¬ 
magnetic  backseat tered  field  [10,11]. 
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In  this  chapter,  the  space-time  integral  equation,  and  how  it 
is  used  to  obtain  a  first  order  correction  to  physical  optics,  is 
discussed  so  that  later,  in  Chapter  III,  its  possible  extension  to 
obtain  higher  order  correction  terms  can  be  presented  clearly.  The 
relationship  between  the  principal  curvatures  and  the  general  [S] 
matrix  is  also  developed  in  Chapter  III.  A  discussion  of  numerical 
verification  with  theoretical  as  well  as  measured  data  is  given  in 
Chapter  V. 


2.2  The  Space-Time  Integral  equation 

An  electromagnetic  wave  incident  on  a  perfectly  conducting 
body  induces  currents  on  the  surface  of  the  scatterer,  which  in  turn 
radiate  and  produce  the  scattered  field.  The  current  distribution 
produces  a  vector  potential  given  by 


*p(r,t)  = 


If 


j(r]x) 


dS ' 


where  J(r,t)  is  the  induced  surface  current  density  at  time  t,  r  is 
the  position  vector  to  the  observation  point,  r1  is  that  to  an  Integra- 
tion  point,  R  =  |r  -  r'|,  t  =  t  -  R/c.  The  geometry  is  illustrated  in 
Fi gure  2. 1 , 

The  total  magnetic  field  H  is  equal  to  the  sum  of  the  incident 


H(r,t)  =  H j  (r,t)  +  H s  ( r , t ) 


where  Hs  (r ,t)  =  V  x 
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=  7j-~  /-■' 's  IiJ(r',r))  X  aR  dS  1  [  1 6 J 

The  operator  L  is  defined  as  1/R^  +  (1/Rc)3/3t,  and  aR  is  the  unit 
vector  of  R.  An  expression  for  J  can  be  obtained  by  shrinking  the 
observation  point  to  a  point  on  the  scatterer's  surface  under  a 
limiting  procedure.  For  the  case  of  a  perfectly  conducting  scatterer, 
the  H-field  boundary  condition  (J  =  an  x  H) ,  with  an  being  the  unit 
outward  normal  vector  to  the  surface,  can  then  be  invoked  to  yield 
3'  vector  integral  equation  as  shown  in  Bennett  [23] 

J(r,t)  =  2an  x  h.  (r ,  t)  +  fJ's  an  x  {IJ(r',T)  x  aR  1  dS1  (2.1) 

The  first  term  on  the  right-hand  side  of  (2.1)  is  identified  a,  the 
physical  optics  approximat ion  and  is  also  the  source  term,  whereas 
the  second  integral  represents  the  contribution  of  retarded  currents 
at  points  on  the  scatterer's  surface  other  than  the  observation  point 
P(r)  i.e. 

J  (r, t)  =  JpQ(r , t)  +  J£  (r , t) , 

with  J£  (r ,  t)  =  anx{LJ(r',x)xaR}dS'  (2.2) 

being  the  correction  to  the  physical  optics  current  3  .  For  the  far 

scattered  magnetic  field  (Hs),  with  R  -*■  r,  it  can  be  shown  that  [8] 

roHs<r,t)  =  ffs  hpf(r'  ,t)]  x  ap  dS '  (2.3) 

where  rQ  is  the  radar  range.  Using  (2. 3).  the  physical  optics  approxi¬ 
mation  for  the  far  scattered  impulse  response  field  was  derived  [8,23] 
and  Is  equivalent  to  the  Kennaugh-Cosgr i f f  formula  [10-12]  which  can 


be  wr i tten  as 


yields  the  first  order  correction  current 
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-T7ffSE[L JU(;,-r)][(an 
+  2V  //s  [iJv(r‘  ,t)][  (an 


aR>au 


aR)av 


(an 

(a„ 


au)aR]  dS 
a^)aR]  dS 


where  the  primed  quantities  are  associated  with  the  integration  point 

and  au,  ay  are  the  unit  vectors  tangent  to  the  principal  curves  at 

the  point  of  interest.  The  geometry  information  embedded  in  a  .  aD 

n  k 

can  be  extracted  by  expressing  aR  =  (r  -  r')/R  in  a  Taylor  series 
expansion  as  follows: 

-*■ 

r1  -  r  =  r  Au  +  r,,Av 
u  v 


P* (r'l 


+  IT  ^uu(Au)2  +  2ruv(Au)(Av)  +  rw(Av)2] 

+  JT  ['uuu(Au)3  +  3ruuv(Au)2(Av)  +  3ruvv(Au) (Av)2 

+  ;vvv(Av)3] 

+  .  .  .  , 

-+■ 

-*■  d  r 

where  etc.  The  above  series  describes  any  neighboring  point 

(Au,Av)  in  the  vicinity  of  r  in  terms  of  the  derivatives  of  r  at  the 

specular  point.  The  geometry  of  the  small  patch  is  thus  extrapolated 

in  terms  of  the  properties  of  the  surface  at  the  specular  point. 

The  scalar  product  of  a„  and  aD  expressed  by  a  series  truncated 

n  tt 

at  second  order  (i.e.  terms  (Au)2  and  (Av)2,  etc.)  introduces  (E,  F, 
G)  and  (L,  M,  N) ,  the  coefficients  of  the  first  and  second  fundamen¬ 
tal  forms  of  the  surface  r(u,v)  [13, 1^].  To  simplify  the  algebra, 
the  principal  curves  are  chosen  as  the  parametric  curves  to  represent 
the  curvilinear  mapping  l*(u,v)  for  the  surface  of  the  patch,  thus 


forcing  F  and  M  to  zero  [8,73,14].  The  principal  curvatures  a  long 
ay  and  ay  at  the  specular  point  are  obtained  as  [13] 
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respectively.  Finally,  LJ  (r‘ ,t)  and  LJ  (r',x)  can  be  approximated 
as  Ju(r,t)/R2  and  Jv(r,t)/R2,  respectively,  by  assuming  that  the 
currents  are  spatially  constant  on  a  small,  flat,  circular  patch 
of  radius  eQ.  With  the  procedure  outlined  above  the  analytic  expression 
arr i ved  at  in  [8]  is 

K, ,  -  Kw 

J£(r,t)  =  [ a u J u ( r , t)  -  avJv(r,t)]  - jj -  Eq  (2.5) 


The  corresponding  first  order  far-field  impulse  response  correction 
was  obtained  by  assuming  physical  optics  currents  for  Ju  and  J in 
the  above  equation  and  then  substituting  it  into  (2.3).  A  crucial 

assumption  made  in  [8]  is  that  the  patch  radius  E  increases  linearly 

o 

with  time  t,  spreading  from  the  specular  point  at  the  leading  edge. 
The  final  expression  for  the  first  order  correction  to  the  scattered 
far  field  is 


Ku  -  Kv 


r  H  ,  ,  (r ,  t)  *  — i — 

o  s(pol)  H  7T 


.  [a  H  .  -  a  H  .]  22. 

u  ui  v  vi  dt 


where  H  .  and  H  .  are  the  components  of  H.  in  the  directions  of  a 
u  i  v  i  i  u 

and  av,  respectively. 

It  is  clear  that  the  first  order  correction  exhibits 
depolarization  effects,  which  are  proportional  to  the  difference 
in  the  principal  curvatures  at  the  specular  point.  Moreover,  the 
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first  order  correction  takes  the  functional  form  of  the  first 
derivative  of  the  silhouette  area  function  A(t),  whereas  the 
physical  optics  far  field,  which  exhibits  no  depolarization 
effects,  takes  the  functional  form  of  the  second  derivative 
of  A(t).  The  practical  aspects  of  using  this  first  order 
correction  in  geometry  extraction  (curvature  difference  at  th 
specular  point  in  this  case)  are  analyzed  in  the  next  chapter 
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CHAPTER  I  1  I 

CURVATURE  RECOVERY  FROM  HIGH  FREQUENCY  [S]  MATRIX  ELEMENTS 

3.1  Derivation  of  the  Phase-Curvature  Relationship 

The  polarimetric  scattering  data  measured  with  a  monostatic 
radar  system  are  given  by 


s  .  s.„ 

'  |S,  |ej*H 

|SlJej^2 

11  12 

1  11  1 

1  121 

S„  s„ 

|S,,  |eJ>21 

|  I eJ '^22 

21  22 

1  21  1 

— 

’  22  1 

The  general  polarization  geometry  with  respect  to  the  principal 
directions  at  the  specular  point  is  shown  in  Figure  3-1-  The  elements 
and  S^2  represent  the  backscattered  signals  when  the  transmitted 
and  the  received  polarizations  are  identical,  i.e.  a1  and  a2» 
respectively.  On  the  other  hand,  and  S^1  represent  the  cases 
when  the  transmitted  and  the  received  polarizations  are  orthogonal 
to  each  other  (transmi  t-a^ ,  receive-a2  S^,  etc.).  In  order  to 
relate  the  measurable  [S]  matrix  to  the  theory  presented  in  Chapter 
II,  the  total  physical  optics  scattered  far  field  (i.e.  the  physical 
optics  and  first  order  correction)  is  transformed  from  the  time  (t) 
domain  to  the  frequency  (k,  the  wave  number)  domain  by  using 
Fourier  transformation.  In  the  time  domain,  combining  (2. A)  and  (2.6), 
the  total  high  frequency  scattered  far  field  is 


Ku  -  Kv  (j 


r  H  (r, t)  =  A ( t )  a„  +  - -  .  y-  A(t) 

os  2  7T  (jj.2  Hi  A  rr  dt 


[(aH.-  au}  au  '  (aH.-  av}  av] 
1  1 


(3.1) 


The  Fourier  transform  of  (3-1).  with  the  initial  condition 


Ap (0)  =  <F.T.[A(t) ])k=0  =  0, 


I i nes  of  principal 
curvatures 


a^,  a^  transmitting  or 
receiving  basis 
vectors 


ure  3.1  Specular  Point  Coordinate  Systems 


r 


yields 
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r<A(7’k)  =  FT  Ok)2AF(k)  aH. 


+  (jk)Ap(k)  [(aH..  au)au 


(aH_.  av5av]  (3-2) 


where  Ap(k)  =  F.T.[A(t)]. 


From  the  geometry  in  Figure  3-1,  it  is  seen  that 


a,  *  sin  a  a  +  cos  a  a  and 
1  u  v 


a„  =  cos  a  a  -  sin  a  a 
2  u  v 


(3.3) 


For  measuring  Sf,  the  transmitter  polarization  becomes  a  =  a,,  and 

I  I  H  |  I 

the  receiver  polarization  becomes  a  =  a..  Thus,  using  (3-2),  one 

n  r  I 

has 


'11 


=  2V  (j  k)  2Ap  (k)  (a  j .  a  ^ ) 


+  (jk)AF(k )J^Uiv  aj  Cay  i0) 

-  (a , .  av)  (a,,  a^) ] 


(3.4) 


Similarly,  for  S22,  aH<  *  a2>  aHr  =  a2>-  for  S^,  aH  *  aJt  aH(.  =  a2; 
and  for  S12,  aH_  *  a 2>  aHr  -  Now  using  (3-2),  (3.3)  and  (3-4), 
one  gets  (ignoring  scale  factors) 

K  *  K 

S11  =  IT  0k)2AF(k)  -  (jk)AF(k)  COS  2a 


K  -  K 


S22  =  IT  0'k)^AF(k)  +  (jk)AF(k)  cos  2a 

K  -  K 

S21  =  (jk)Ap(k)  -~jj-  ---  sin  2a 


(3.5) 

(3-6) 

(3.7) 


^21  =  ^12  (reciprocity  satisfied) 
Using  (3-5)  and  ( 3 . 6 ) ,  it  can  be  shown  that 


K  -  K 
u  v 


S  -  S 

k  11  a22 

cos  2a  Sj 1  +  S22 


cos  2a 


1  -  j? 
1 


(3.8) 


where  $  =  /?eJ'^8 


^4eJ(*22  '  *11} 

I  i  I 


It  is  clear  that  in  order  for  the  equality  in  (3.8)  to  hold  true 
and  therefore  represent  a  physical  case  in  which  («u  -  K  )  is  a 
real  number,  one  needs 


i  _  Ji  j 

Re  - —  =0,  implying  1  -  /r  -  0 

1  +  R 


where  Re  stands  for  real  part. 

Thus  for  (3-8)  to  represent  a  physical  situation,  the  condition 
required  is 


TsTT 


=  1,  implying  |S  |  =  |S 


From  electromagnetic  scattering  theory,  it  is  known  that  the  above 
condition  is  attained  at  relatively  high  frequencies  (i.e.  physical 
optics  to  geometrical  optics  region).  It  is  interesting  to  note 
that  algebraic  manipulation  of  (3-8)  independently  points  out  that 
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it  is  a  high  frequency  formula  which  is,  of  course,  true  since 
the  physical  optics  approximations  for  the  scattered  fields  are 
being  used. 

With  the  condition  (3.9),  now  (3-8)  can  be  written  as 


K 

u 


2 


K 

v 


k 

cos  2a 


tan  (4>d/2) 


Where  4>d  =  i>22  *  <)1 ,  • 


(3-10) 


In  the  rest  of  the  text,  the  above  expression  will  be  referred 
to  as  the  "phase-curvature"  relationship. 

For  the  inverse  scattering  applications  a  i s  an  unknown 
quantity  (Figure  3.1).  Thus,  before  (3.10)  can  be  used  to  recover 
the  curvature  difference  at  the  specular  point,  a  needs  to  be 
determined.  For  this  purpose,  consider 


from 


K 

u 


(jk)  - 


(3.5)  and  (3.7) 


_  K 

v  .  „ 

-  s  i  n  2a 

_  __ 

u  v  « 

- - -  cos  2a 

.  Using  ( 3 - 1 0 ) ,  one  obtains 


S21  (tan  2a)  (tan  (i^/2)) 

51 1  tan  (<f>  ,/2)  +  j 

d 

-S 

a  =  i  tan  1  {  ~ —  (1  +  j  cot  U  ,/2))  }  (3-11) 

L  bn 


Once  again  for  a  to  be  a  real  angle  (and  therefore  representing  the 
given  physical  situation)  one  needs 


D  (1  -  j  cot  (<f,/2)) 
o  d 
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where  is  a  real  number. 

Using  (3.1l),  without  applying  any  unitary  transformation, 
the  cross-polarized  nulls  of  a  given  scattering  matrix  are  known. 

Two  special  cases  of  (3-10)  are  when  a  =  0,  or  tt/2  (i.e.  the 
incident  linear  polarization  coincides  with  one  of  the  directions  of 
the  principal  curvatures  at  the  specular  point).  For  these  special 
cases,  there  is  no  depolarization  of  the  energy  in  the  backscattered 
direction  =  0).  The  corresponding  [S]  matrix  is  referred 

to  as  the  cross-pol .  null  scattering  matrix.  The  phase-curvature 
relationship  becomes 


K  -  K 
u  v 


2 


+  k  tan 


(3- 12) 


In  chapter  V,  the  validity  of  (3-12)  will  be  numerically  tested 
with  both  theoretical  and  experimental  scattering  matrix  data. 

3-2  Numerical  Analysis 

In  this  section,  a  numerical  analysis  of  applying  the 
phase-curvature  relationship  on  the  poiarimetric  scattering  matrix 
of  a  prolate  spheroidal  scatterer  is  presented.  The  2:1  prolate 
spheroid  was  used  as  a  test  case  because  of  its  well-defined  finite 
curvature  difference  at  any  point  on  the  surface  and  also  because 
the  theoretical  as  well  as  experimental  data  over  a  large  range  of 
frequencies  were  readily  available  for  this  object.  For  the  time 
being,  since  both  theoretical  and  experimental  data  are  available 
only  for  the  special  cases  in  which  the  incident  polarization 
coincides  with  one  of  the  directions  of  the  principal  curvatures 
at  the  specular  points  (i.e.  points  on  the  equator  of  the  prolate 
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spheroid  (Figure  3-2)),  equation  (3-12)  rather  than  the  more  general 
phase-curvature  relationship  (3-10)  is  directly  tested.  In  these 
cases,  the  incident  polarizations  are  along  the  vertical  and  horizontal 
directions.  To  verify  the  theory  of  curvature  recovery,  it  is  essential 
to  check 

(i)  whether  the  right-hand  side  of  (3-12)  will  approach  the  actual 
value  of  (K^  -  K  )/2  for  a  given  equatorial  specular  point  as 
k  increases; 

(ii)  whether  k  Re  {  (1  -  ??)/(!  +  r?)  }  will  tend  to  zero  despite  that 
k  increases; 

(iii)  whether  the  imaginary  part  of  k  {  (1  -  $)/(l  +  R)  }  will 

settle  to  the  constant  value  of  (K  -  K  )/2; 

u  v 

A  more  interesting  and  compact  presentation  of  results  is  to 

plot  the  right-hand  side  of  (3.8)  on  a  complex  plane  (i.e.  the 

imaginary  part  versus  the  real  part  of  f  k  ( 1  -  $) / ( 1  +  R)  }.  It 

is  predicted  that  this  scattering  chart  will  be  a  spiral  which,  as 

frequency  increases,  will  converge  to  (or  hover  around)  a  point  on 

the  imaginary  axis.  The  distance  of  this  point  on  the  imaginary 

axis  from  the  origin  will  be  equal  to  the  required  value  of  (K  - 

K  )/2  for  the  specular  point  of  interest.  It  is  expected  that  this 

complex  plane  plot  will  be  particularly  useful  when  the  input  data 

are  not  very  accurate  and  are  noisy  (measurement  data).  In  the  above 

tests,  the  exact  value  of  the  curvature  difference  K  -  K  has  been 

u  v 

calculated  by  using  Minkowski's  support  function  for  ellipsoidal 
surfaces  [17].  Another  way  is  through  the  use  of  differential 
geometry  (Appendix  l).  The  value  of  k  is  normalized  with  respect 
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to  the  length  of  the  semi-minor  axis  of  the  prolate  spheroid. 

It  is  useful  to  consider  the  following  types  of  measurement 
errors  for  the  interpretat ion  of  the  numerical  results  and  the 
explanation  of  the  deviation  from  theoretical  prediction  given  in 
chapter  V. 


3.2.1  Relative  Phase  Error  between  TE  and  TM  Incidences 

If  the  target  is  not  illuminated  simultaneously  with  TE 

(vertical)  and  TM  (horizontal)  polarizations,  then  even  an  offset 

of  one  millimeter  between  the  separate  positions  of  the  same  antenna 

along  the  direction  of  incidence  will  cause  phase  distortion  of 

about  ten  degrees  to  the  phase  difference  (<i>vv  ”  ®^)  at  a  frequency  of 

A  GHz.  If,  however,  the  target  is  slightly  displaced  in  the  direction  of 

propagation  with  the  TE  and  TM  measurements  simultaneously  made,  then  both 

t>  and  i>.,  are  distorted  to  the  same  extent.  Thus  the  relative  phase 
vv  hh 

difference  and  the  subsequent  calculations  are  not  affected  at  all. 

To  investigate  how  the  scattering  chart  is  affected  assuming 
that  a  relative  phase  error  of  e  radians  does  occur,  the  complex 
ratio  (l  -  ??)/(l  +  R)  can  be  broken  down  into  its  real  and  imaginary 
parts,  with  <f>^  replaced  by ^  +  e) . 


1 

+  * 
1  -  2 

m  - — 


De 

-2R  sin  (  4>  .  +  e) 

_ d _ 

De 


where  De  =  1  +  R  +  1R  cos  (i>,  +  e) .  Note  that  a  spheroidal  wave  function 

d 

expansion  could  be  used  for  error  estimation. 

It  is  obvious  that  the  real  part  is  much  more  resistant  to 


changes  caused  by  the  error  e  than  the  imaginary  part.  Thus  the 


spiral  will  become  a  helix  elongating  mainly  along  the  imaginary 

axis  and  away  from  the  exact  value  of  (K  -  K  )/2.  In  Chapter  V, 

u  v 

the  relative  phase  error  will  be  simulated  within  the  set  of 
theoretical  data. 

3-2.2  Rotation  of  Target  with  respect  to  Incidence  Direction 
If  the  target  is  rotated  with  respect  to  the  direction  of 
incidence  (i.e.  the  incident  polarizations  are  not  along  the 
directions  of  principal  curvatures),  then  it  is  the  general 
phase-curvature  relationship  (3-10)  which  should  be  tested  instead. 
Hence,  a  mu  1 1 i p I ica t i ve  factor  of  cos  2a  (Figure  3-l)  will  account 
for  the  sole  effect  of  rotation  error.  The  scattering  chart  will 
retain  its  spiral  shape  however. 

3.2.3  Canting 

If  the  target  is  slightly  canted  with  the  z-axis  still  being 
horizontal,  then  obviously  the  effect  will  be  that  which  resulted  from 
changing  the  aspect  angle  $  (Figure  3-2). 

3.3  Discussion  of  Second  Order  Corrections  to  the  Physical 
Optics  Approximation 

The  first  order  correction  to  the  physical  optics  approximation 

due  to  Bennett  et  al  [8]  has  been  shown  to  be  proportional  to  the 

difference  in  principal  curvatures  and  to  have  the  functional 

form  of  the  first  derivative  of  the  silhouette  area  of  the  target. 

It  might  be  expected  that  more  geometry  of  the  specular  point  can 

be  identified  by  extending  the  first  order  correction.  For  a  general 

case  in  which  (K  -  K  )  is  not  zero,  Geometrical  parameters,  such  as 
u  v 
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derivatives  or  mixed  derivatives  of  principal  curvatures,  torsion 
of  the  principal  curves,  etc.,  can  appear  in  the  second  order 
correction  term.  Another  motivation  in  the  pursuit  of  the  second 
order  correction  is  the  possible  relaxing  of  the  high  frequency 
restriction  required  for  the  validity  of  the  first  order  correction 

term.  In  order  to  elaborate  on  this  conjecture,  note  that,  since 

d  .  9 

— -  of  the  physical  optics  term  corresponds  to  (j k) ^A. (k)  in  the 
dt^-  r 

dA 

frequency  domain,  and  of  the  first  order  correction  term 
corresponds  to  (jk)A[_(k),  there  might  be  a  trend  that  higher  order 
corrections  become  more  dominant  (important)  for  lower  frequencies 
within  the  high  frequency  realm  of  physical  optics. 

One  approach  to  the  problem  of  extension  is  to  take  into 
account  the  higher  order  derivatives  of  r(u,v)  that  were  truncated 
in  the  first  order  correction.  The  first  order  term  has  been 
obtained  by  retaining  the  second  order  derivatives  of  r  with  respect 
to  u  and  v.  To  obtain  the  second  order  correction  term,  the  integrand 
in  (2.2)  for  J£  has  been  expanded  with  inclusion  of  the  higher  order 

0  -V  # 

surface  derivatives  of  r.  After  some  algebraic  manipulations,  the 

integrand  was  written  as  a  sum  of  several  terms  like  a  ,  a  ,  r  , 

3  u  v  uu 

r  ,  ...  (up  to  the  third  order  derivatives)  multiplied  by  the  powers 
of  Au  and/or  Av.  This  is  expected  since  these  vectors  introduce  the 
higher  order  depolarization  concepts.  In  contrast,  the  first  order 


correction  term  is  in  the  directions  of  a  and  a  explicitly. 

u  v 

—V  — ►  "  ' 

The  vectors  such  as  r  ,  r  ,  etc.  can  be  resolved  in  a  ,  a 

uu  uv  u  v 

and  an  directions  by  introducing  the  Christoffel  symbols  through  the 
use  of  the  Gauss  and  Weingarten  equations  in  differential  geometry 
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[12,13]-  These  Christoffel  symbols  depend  on  the  basic  geometric 
parameters  of  the  surface,  namely,  (E,  F,  G)  and  (L,  M,  N) .  All 
terms  in  the  a  direction  can  be  neqlected  on  the  basis  of 
physical  considerations  that  the  induced  surface  current  can not 
have  a  component  in  the  normal  direction  to  the  surface.  The 
appearance  of  the  a^  term  is  due  to  the  truncation  of  the  Taylor  series. 
The  totality  of  all  terms  in  the  series  should,  in  theory,  nullify 
the  current  component  in  the  a^  direction.  Since  the  depolarization 
derivative  vectors,  and  hence  the  Christoffel  symbols,  are  all 
evaluated  at  the  specular  point,  they  can  be  taken  outside  the 
surface  integral  in  (2.2),  which  then  is  written  as  a  sum  of 
integrals  of  the  type 

//.  -^!_dS'  ,  //  iM_  dS'  ....  etc. 

R-5  rh 


For  direct  scattering  problems  the  integration  may  yield  complicated 

analytic  solutions  despite  the  fact  that  S£  is  known.  For  inverse 

scattering  problems,  where  Se  is  unknown,  the  integration  may  be 

approximated  by  assuming  a  flat  patch  for  S£.  Since  the  Taylor 

series  dictated  a  small  patch,  leading  edge  effects  and  thus  the 

high  frequency  restrictions  must  be  adhered  to.  An  approximate 
■+ 

expression  for  J£  can  thus  be  obtained  by  choosing  a  circular 

patch  of  very  small  radius  e  as  the  most  simple  case.  However, 

o 

the  expression  thus  obtained  has  been  found  to  have  a  factor  of 

e  ^  in  contrast  to  e  in  the  first  order  correction.  The  relatively 

o  o 

small  value  of  eJ  compared  with  e  renders  the  expression  for  the 

o  o 

second  order  correction  current  thus  obtained  insignificant.  The 
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fa r- field  correction  term  corresponding  to  this  second  order 
current  term  then  has  very  little  significance  in  comparision 
to  the  physical  optics  or  the  first  order  term,  and  therefore 
is  of  little  consequence  to  most  practical  situations. 

Nevertheless,  it  is  suggested  that  one  may  assume  some 
simple  known  curved  patch  (instead  of  a  flat  patch),  and  take 
into  account  the  whole  vicinity  of  the  specular  point  rather 
than  only  the  point  itself,  in  further  pursuit  of  specular 
geometry  through  the  backscattered  leading  edge  returns. 
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CHAPTER  IV 

THE  HF  CURVATURE  RECOVERY  MODEL  AND  THE  TRANSFORMATION  INVARIANTS 
OF  THE  SCATTERING  MATRIX 

A. 1  Relation  to  the  Transformation  Invariants  of  [S] 

The  scattering  matrix  [ S  (A , B ) ]  with  respect  to  an  orthogonal 
basis  (A,B)  can  be  transformed  to  tS’(A',B')]  with  respect  to  another 
orthogonal  basis  (A 1 , B 1 )  through  a  unitary  transformation  [18,19. 
20],  The  transformation  is  invariant  and  satisfies  [3,38] 

Span  { [S (A, 8) ] }  *  Span  {[S'  (A ' , B ‘ ) ] }  =  invariant  (4.1) 

and  Det  {[S(A,B)]}  =  Det  { [ S ' (A1 , B 1 ) ] }  =  invariant  (A. 2) 

where  Det  stands  for  the  determinant  of  the  scattering  matrix  and  the 
span  is  defined  in  the  following. 

Applying  equations  (3-5  ~  3-7)  for  a  basis  (A,B)  with 
polarization  angle  a  (Figure  3-1), 

Span  { [S (A, B) ] }  =  | SAA | 2  +  I $BB i 2  +  2 | SAB ) 2 

=  (~)2k2lAF(k)  1  2 1  j  k  -  ~  2-  -V  cos  2a  |  2 

K  -  K 

+  (~)  2k2  |Af  (k)  |  2[  jk  +  —  2  V  cos  2a|2 
K  k 

+  2(^-)2k2lAp(k)  [2[— y-V  sin  2a]2 
=  2(^)2k2|AF(k)|2{k2  +  [JL-J^-]2}  ,  (A. 3) 

which  is  independent  of  a. 

For  a  given  frequency  (high  enough  so  that  the  theory  is 
valid)  and  a  given  aspect,  the  right-hand  side  of  equation  (A. 3) 


is  indeed  invariant. 
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Similarly, 

f[S<A,B)]}.S4ASeB-SAB2 

i  ,  ,  ,  ,  K  -  K  2  . 

=  (jjf)  (jk)  Ap  (k)  {  ( j  k)  -  L-^— 2 — — 1  (cos  a 

OO  19  99  K  -  K  2  - 

-  sin  a)  }  -  (— )  (jk)  Ap(k)[-^— 2 — -]  s  i  n  2a 

=  (^7)  k2A‘(k){k2  +  [^-2  V]  }  (4.4) 

which  is  also  invariant  and  independent  of  a. 

Hence  equations  (4-3)  and  (*4.4)  are  the  high  frequency  versions 
of  the  invariance  equations  (4.1)  and  (A. 2). 


4.2  The  HF  Scattering  Ratio 

It  is  interesting  to  define  D  (referred  to  as  "the  HF 
scattering  ratio"  in  this  thesis)  as  the  ratio  of  Det  {[S(A,B)]} 
to  Span  { [S (A, B) ] } , 


D  = 


S  S 
AA  BB 


AB 


ISAa!2  +  iSBBi2  +  2'SABI 

Ap2(k) 


2 | Ap (k) \‘ 
j  20 


=  0.5  e 


(4.5) 


where  0  =  Arg  A^(k)  =  the  phase  of  Ap(k)  (4.6) 

For  the  case  in  which  incidence  polarization  is  along  one 
of  the  directions  of  principal  curvatures,  |S  |  equals  zero.  Let  S^A 
|SAA|e^AA  and  so  on,  thus 


be 
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D  = 


lSAA‘ I SBB i  e 


j  (<fAA  +  ^Bb) 


iSAA^  +  lSBB^ 


Thus  6  =  Arg  Ap(k)  = 


^AA  +  ^BB 


(4.7) 


and 


isaa!  IsbbI  =  j_ 

I-  | 2  | _  i2  ~  2  ,  which  requires  that  | S  j  =  jS  j.  This 

|5AA'  +  ' SBB '  AA  BB 


condition  is  consistent  with  high  frequency  electromagnetic  scattering. 
Returning  again  to  the  general  case  in  which  |S^_|  4  0, 


S  S  -S  2  =  I  s  Ms  leJ'  ^AA  + 
AA  BB  AB  1  AA11  BB 1 


is  i2J2^ab 

1  ^AB 1  e 


=  ( i SAA I  IsbbIcos  (*AA  +  4>bb)  -  i sab i 2cos  2W 


AB1 


+j  ( i SAA I  I SRB I  5 1 n  (^AA  +  *BB)  ‘  I  SAB  |2sin  2d>/x3) 


AA1 1  BB 


AB1 


Considering  amplitude  only, 


{  |  S  |  2  I  S  I  2 

1  AA1  BB1 


2IS..I  |S„[  |S..|2  cos  (iAAHB6-2*AB)  *  Is.j'1  )J 


AA 1 1  BB 1 1  AB1 


AB1 


AA  1 


+  ISBb'Z  *  2 1 SAB 1 2 


The  above  equation  is  an  identity,  if 


(i) 


AA1 


=  S 


BB1 


(4.8) 


(ii) 


AA 


BB 


-  2* 


AB 


(4.9) 


simultaneously  hold. 

Again,  (4.8)  is  consistent  with  high  frequency  scattering. 
It  should  be  noted  that  when  (jk)  is  neglected  compared  to 
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2 

( j  k )  in  (3-5)  and  (3-6)  for  high  frequencies,  both  (4.  8)  and 
particularly  (4. 9)  result.  Hence,  in  general, 

| 0 |  =  0.5  (4.10) 

is  not  a  trivial  result  merely  from  (4.8)  as  in  the  previous 
special  case,  but  rather  a  consequence  of  the  first  order  correction 
to  physical  optics. 

4.3  Interpretation  of  the  HF  Phase  Sum  (<j>^  +  <j>  ) 

From  the  above  analysis,  with  only  the  phase  being  considered, 
it  is  found  that 


tan  29  =  tan  2(Arg  A^(k)) 


|saaMsbbI  s!"  4AA 


BB 


)  - 


SAB'  S!n  2*AB 


ISAaI iSBBl  COS  (*AA  +  V  " 


IsabI  cos  2Vb 


Assuming  (^.S)  and  neglecting  js^l  compared  to  |$AA|,  the  above 

equation  becomes  (4-7),  which  is  now  also  valid  for  the  general 

case.  Its  validity  enables  the  phase  sum  (4>An  +  <t>DD)  to  be  inter- 
j  A  A  dd 

preted  as  twice  the  argument  of  the  Fourier  transform  of  the 
s i Ihouettei area  of  the  target  within  the  high  frequency  range. 


4.4  Numerical  Analysis 

If  the  real  and  imaginary  parts  of  D,  the  scattering  ratio, 
are  plotted  on  a  complex  plane,  a  circle  of  radius  0.5  will  be 
expected  for  high  frequency  polarimetric  data  input.  Since  phase 
changes  rapidly  with  frequency,  a  circle  rather  than  a  cluster  of 
points  of  phases  about  2  Arg  Ap(k)  will  appear.  This,  and  a  direct 


plot  of  | D |  versus  k  will  be  shown  in  Chapter  V.  One  significance  of 
these  plots  is  to  check  the  accuracy  of  experimental  polarimetric  data 
of  high  frequencies  by  observing  the  deviation  from  the  perimeter  of 
the  circle  of  radius  0.5. 

4.5  Transformat  ion  to  Circular  Polarization  Basis  Vectors 

The  orthonormal  vectors  along  the  horizontal  and  vertical 
directions  are  usually  chosen  as  the  polarization  bases  (denoted 
by  (H,V))  for  both  the  transmitting  and  the  receiving  systems.  However, 
a  circular  polarization  basis  pair  may  also  be  used  [36,40],  particularly 
in  radar  meteorology,  in  which  circular  polarization  has  a  particular 
appropriateness  on  account  of  the  direct  correspondence  between  the 
mean  orientation  angle  and  the  relative  phase  of  received  circular 
polarization  components [ 29] .  Circular  polarization  has  also  been  utilized 
in  the  backscatter  measurements  of  dielectric  spheroids  [32,33].  One 
way  to  investigate  the  form  which  the  phase-curvature  relationship  may 
take  in  circular  polarization  basis  is  to  transform  [S(HV)]  with  respect 
to  the  linear  basis  (H,V)  to  [C(RL)]  with  respect  to  a  circular  polarization 
basis  (R.L)  [31,6,27].  Such  transformation  of  [S]  depends  on  the 
specification  of  the  transformation  of  (H,V)  to  (R,L)  by  a  matrix  [ T ( RL ; HV ) ] 

R  H 

=  [ T ( RL ; HV ) ]  (4.11) 

Lj  V 

where  L  and  R  denote  the  left-circular  and  right-circular  polarization 
vectors,  respectively.  The  left-circular  and  right-circular  senses  are 
defined  in  Figure  4.1.  tn  general,  the  transformation  of  the  linear  basis 
(H,V)  to  any  other  orthonormal  basis  (A,B)  (not  necessarily  circular 
polarization  basis)  through  [ T ]  must  satisfy  the  normalization  requirements 
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H  .  H  =  1 

J. 

V  .  V  =  1 

J. 

A  .  A  =  1 

B  .  B  =  1 


and  the  orthogonality  requirement 

.t. 

H  .  v  =0 

.U 

A  •  B  =0 


These  requirements  can  be  shown  to  be  mathematically  equivalent  to 


[T]*T  =  [T] 


which  satisfies  the  definition  of  a  unitary  matrix.  Hence  [T]  is  a  unitary 
matrix,  and  its  most  general  form  can  be  written  as  [30] 


[T] 


eJ'^1  cos  B  e^2  sin  (s 

-e^3  sin  B  e^**  cos  B 


with  (*>2  "  ^1  =  ^4  ”  l*>3'  The  most  general  basis  (A,B)  is  an  elliptic 
one.  When  all  phases  are  set  to  zero,  [T]  is  just  an  ordinary  rotational 
matrix  which  rotates  (H,V)  to  another  linear  basis  by  an  angle  6.  An 
example  is  given  by  the  invariant  transformation  described  by  equations 
(4.3)  and  (4.4),  which  shows  from  the  curvature  recovery  model  that 
rotational  transformat  ion  alone  renders  the  invariants  independent  of 
the  polarization  angle  (Figure  3-1).  A  more  general  case  can  be  given 


by,  for  instance,  setting 


D  =  n/4 


=  ^  =  0 
4>2  =  ~4>3  =  tt/2 


The  corresponding  [ T ]  then  becomes 


[T]  = 


72 


(4.12) 


and  the  corresponding  (A,B)  reduces  to  a  circular  polarization  basis 


R 

"l  j  ‘ 

H 

1 

"  72 

L 

j  1  _ 

V 

By  limiting  the  transformation  from  (H,V)  to  circular  polarization 
basis  only,  any  polarization  vector  can  then  be  expressed  in  terms 
of  either  basis.  For  instance,  the  incident  electric  field  polarization 
vector  E'  can  be  written  as 


E' 


EH 


H 

V 


R 

L 


where  the  linear  phasor  components 


can  be  related  to  the 
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circular  phasor  components 


E 


i 

R’ 


EL3S 


f o  t ! ows 


(4.13) 

where  the  symbols  *  and  [  denote  conjugation  and  t ransDos i t i on, 
respectively.  In  (4.13),  the  incident  polarization  vector  can  be 
regarded  as  being  transformed  in  changing  the  polarization  basis  as  speci¬ 
fied  by  the  unitary  matrix  [T]  given  by  (4.11).  The  scattered 
polarization  vector  £s  can  similarly  be  transformed.  E_S  and  E_'  can 
both  be  specified  in  terms  of  the  circular  polarization  basis.  Vet, 
it  is  preferable  to  use  distinct  systems  of  unit  vectors  to  specify 
the  incident  and  scattered  fields,  so  that  right-hand  elliptical 
polarization  may  have  the  same  sense  with  regard  to  the  coordinate 
system  for  incident  radiation  as  it  does  with  regard  to  the  coordinate 
system  for  scattered  radiation  [30, 7]  (Figure  4.1).  If  (4.11)  is  prescribed  for 
the  incident  system,  the  desired  similarity  of  sense  for  the  two 
coordinate  systems  can  be  accomplished  by  writing 


V" 

E,', 

R 

X 

H 

< 

=  [t(rl;HV)]“ 

Ei 

R 

s 

X 

H 

=  [t (RL;HV)  ]"' 

L 

V 

(4.14) 


for  the  scattered  system.  Because  of  the  conjugation  now  introduced, 
the  relative  phases  in  [T]  are  negated,  and  thus  the  sense  of  rotation 
is  reversed.  Also,  it  follows  that 


4! 

i 


L  3 T  .  •  TSt-  ■ 


/ 


/ 


Right-c 
w. r . t . 


rcular  sense 
scattered  system 


Right-ci rcular 
sense  w. r. t. 
incident  system 


/ 


system 


w.r.t. -  with  respect  to 


Figure  4.1  Incident  and  Scattered  Coordinate  Systems 
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On  adopting  this  choice  for  the  circular  polarization  basis  of  the 
scattered  field,  (4. 15)  follows  from  (4.14),  with  the  double  conjugation 
being  ignored. 


[T (RL; HV) ] 


(4.15) 


in  the  same  way  (4.13)  follows  from  (4.11). 

In  terms  of  the  linear  basis  (H,V),  the  scattering  matrix 
[S (HV) ]  fully  describes  the  scattered  depolarized  field,  with  the 
incident  field  given,  by 


Eu 

E1 

H 

=  [s  (HV)  ] 

H 

E* 

E1 

V 

V 

- 

(4.16) 


Similarly,  in  terms  of  the  circular  basis  (R,L),  it  follows  that 


-  - 

esb 

V 

R 

=  [C (RL) j 

R 

E? 

e! 

L 

L 

(4.17) 


Using  the  definition  of  the  unitary  matrix,  it  follows  from 
(4.16),  (4.17),  (4.13)  and  (4.15)  that 

[C (RL) J  =  [T(RL;HV)j  [S(HV)j  fT(RL:HV)]T  (4.18) 


which  transforms  [S]  to  [C]  through  [T]  defined  in  (4.11),  and  which 
is  of  the  form  of  congruence  transformation. 
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*♦.6  Curvature  Recovery  from  the  Circular  Polarization  [C] 

If  the  unitary  matrix  in  (4.12)  is  adopted,  then  (4.18) 
becomes 


[C (RL) ] 


S  -  S 
HH  VV 


+  j  s 


HV 


s  +  s 

SHH  VV 


S  +  S 
HH  VV 


VV  HH 


+  j  S 


HV 


(4.19) 


provided  that  S  =  S  ,  which  is  true  for  the  monostatic  reciprocal  case 
HV  VH 

(also  CDI  =  C._).  If  conjugation  in  (4.14)  is  used  to  preserve 
the  similarity  of  sense  in  both  the  scattered  and  incident  systems, 
then  an  examination  of  (4.18)  reveals  that  reciprocity  is  satisfied, 
i.e.  [C]  is  symmetric(if  [S]  is  symmetric).  A  different  choice  of  [T] 
may  result  in  a  different  [C].  For  instance, 

r,  ♦/ 


if 


[T] 


1 

72 


1 


(4.12') 


then  [C] 


S  -  S 

HH  VV  .  .  P 
- 2 -  +  J  SHV 


s  +  s 

HH  VV 


S  +  S 
HH  VV 


s  -  s 

HH  VV 


j  s 


HV 


(4.19') 


Returning  to  the  [C]  in  (4.19),  in  view  of  (3-5)  to  (3-7), 


the  matrix  elements  can  be  written  as 
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K  -  K 

k  A_(k)  - ■ - —  [sin  2a  +  j  cos  2a] 

F 


(4.20) 


K  -  K 

On  =  -  k  A^(k)  — -- -  [sin  2a  -  j  cos  2a] 


(4.21) 


CLR  '  CRL  *  -  J  Kk2  V“> 


(4.22) 


Hence , 


K  -  K  2 


2  CLL  CRR 


(4.23) 


which  is  an  equation  of  curvature  recovery  from  the  scattering 
matrix  in  circular  polarization. 

2 

Expressing  the  ratio  C^L  CRR  /  in  terms  of  the  linear 
polarization  [S]  elements, 


CLL  CRR 


(i  -  3) 2  ♦ 


(1  + 


(4.24) 


Combining  (4.23)  and  (4.24), 


K  -  K 
u  v 


♦  jk  {fttisB,* 

CLR 


(4.25) 


1(1  ~  ??)2  +  4(J^V 

\5hh  / 


(1  +%)' 


(4.26) 


Comparing  (4.26)  to  (3.8),  it  can  be  observed  that  the  unknown 


4o 


polarization  angle  a  in  (3.8)  is  being  disguised  in  (4.26)  and  appears 
in  the  form  of  which  is  incorporated  into  the  square  root.  Thus 

one  advantage  of  utilizing  circular  polarization  over  utilizing  linear 
polarization  for  curvature  recovery  is  that  the  polarization  angle 
does  not  have  to  be  determined,  but  still  the  entire  scattering  matrix 
has  to  be  measured. 

It  can  be  seen  that  for  the  case  in  which  Sm,  =0  (a  =  0  or 

nv 

90  degrees),  (4.26)  reduces  to  (3-8)  and  then  (3.12).  As  in  the 
derivation  of  the  phase-curvature  relationship,  the  imaginary  part 
of  the  square  root  of  CR^  C^L  /  CRL  should  give  better  curvature 
difference  and  the  real  part  should  vanish,  as  frequency  is  increasing. 

It  is  to  be  noted  that  under  suitable  conditions,  quantities 
such  as  the  radar  cross  sections  oRL>  aRR,  and  quantities  derived 
from  them  (e.g.  (1  -  $)/(l  +  $) )  yield  meaningful  measurables  in 
measurements  of  the  backscatter  of  dielectric  spheroids  and  hydro¬ 
meteors  [32,33,34];  in  radar  target  discrimination  techniques,  the 

-  |C„.  |  *■  in  [36],  in  view  of  (4.20)  to  (4.22),  can 
1 2i-4 1  A 


quantity  |CRR||CLL,  , -RL 


be  interpreted  as  [- ( 1 /2tt ) |Ar (k) j *] ,  which  reveals  area  information 
for  a  smooth,  convex,  conducting  target  at  high  frequencies. 


4.7  Transformation  Invariants  of  the  Scattering  Matrix 

The  transformat  ion  of  [S]  to  [C]  due  to  transforming  the  linear 
basis  (H,V)  to  the  circular  basis  (R,L)  can  be  achieved  by  introducing 
the  appropriate  relative  phase  between  the  two  orthonormal  vectors  of 
the  (H,V)  basis  in  addition  to  rotation.  To  investigate  if  the  transformat  ion 
is  still  invariant  in  changing  linear  to  circular  bases,  the  Det  of  both 
sides  of  (4.18)  is  taken  : 


-1 


Det  { [C] >  =  Det  { [T] }  .  Det 

=  Det  { [ T  ] }  •  Det 


Det  {[Tj 

1) 

|  Det  UT] 

nr 

=  ej2Arg (Det{ [T] }) 


{[S]l  .  Det  { [ T ]  } 

US]}  .  Det  {[T]* 
Det  { [ S  ] } 


Det  { [ S  ]  > 


(4.27) 


From  the  general  form  of  [T]  given  on  page  33, 

Arg  (Det  { [ T ] } )  =  ^  +  <p^  =  4>2  +  ^ 

Hence, 

Det  { [C  ] }  =  ej2(<f>1  +  4>1*)  .  Det  { [  S  ] }  (4.28) 

It  can  be  easily  proven  that  | Det  { [ T ] } |  =  1  for  the  unitary  matrix 
[  T]  ,  but  it  is  generally  not  true  that  the  determinant  of  a  unitary 
matrix  is  purely  real.  Thus,  the  determinant  of  [ C]  is  strictly 
invariant  iff  Det  { [T] }  is  purely  real,  i.e. 

Det  { [T] }  =  +  1  ,  (4.29) 


otherwise  the  invariant  differs  only  by  a  phase  shift  of  2(<>^  +  i>^) . 
For  examples,  the  [T  ]  given  by  (4.12)  renders 


CRR  CIL  '  CRL 


s  s 
HH  VV 


(4.30) 


as  is  also  evident  from  (4.19).  The  invariant  value  is  the  same  as  that 
in  (4.4)  derived  from  the  curvature  recovery  model,  since  the  rotational 
matrix  possesses  a  real,  unity  determinant.  On  the  other  hand,  the  [T  ] 
in  (4. 12')  will  not  preserve  strict  transformation  invariance,  as  is  also 


evident  from  (4.19')  (i.e.  the  determinants  of  [C]  and  [S]  are  equal 
in  magnitude  but  of  opposite  sign). 
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To  show  from  the  curvature  recovery  mode)  the  invariance  of 
Span  { [  C  ]  } ,  equations  (4.  20)  -  (4. 22)  and  (3-  5)~  (3-  7)  can  be  used  to  give 

Span  {[C] }  =  |CRR|2  +  |CLLj2  +  2  |CRL | 2 


=  Span  { [ S ] >  , 


which  equals  the  invariant  value  in  (4.3). 

It  can  be  shown  that  the  span  is  in  general  invariant,  regardless 
of  whether  (4.29)  holds  or  not.  Denoting  the  trace  of  a  square  matrix 
(i.e.  the  sum  of  the  diagonal  elements)  by  Tr,  it  follows  that 

Span  {[  C] }  =  Tr  { [C ] "[C ] } 

=  Tr  {[ [T][S][T]T]*[ [T][S][T]T]}  (from  (4.18)) 

=  Tr  { [  [T  ]*  [S  ]'"  [ S  ]  [T  ]T  ] } 

[S]  [S]  is  identified  as  the  power  scattering  matrix  [P]  of  [S] 
[39,38],  i.e. 

[P]  =  [S]*[S]  (4.31) 


Let  [P‘]  be  the  power  scattering  matrix  of  [C],  i.e. 

[P' ]  =  [C]*[C]  (4.32) 

=  [T]"[S]"[S][T]T  (4.33) 

Hence, 


Span  { [ C  ] } 


Tr  { [P '  ] } 


*3 


»  Z  PI. 

i  i 
i 


■  i  ?  {Pki  I?  Vi,i> 

k  1  i 


Since  the  set  of  conditions 

,2  , 


It  | ^  . 

1  T 

it,,i  + 

•T21 

It,2I2* 

!T22 

JL 

T  .  T ,  _  + 

T 

1112 

21 

is  equivalent  to 

the  - 

some  algebraic  manipulation)  it  follows  that 


Span  { [C  ] }  =  Pn  +  P22 


Thus , 


Span  {[C]  >  =  Tr  { [P '  ]  >  =  Tr  f [P] }  =  Span  { [S] > 


(MM 


fr.8  An  Interpretation  of  the  Scattering  Ratio 

The  radar  cross  section  a  has  been  defined  in  [ ^ 1 ]  as  follows  : 

or(.  =  |hr  •  [SlhV  (4. 35) 

where  h*  is  the  transmitting  polarization  vector  and  _hr  is  the  antenna 
height  [20].  In  this  definition,  h_  and  h_  are  normalized  to  unity.  It 
has  been  shown  in  [20]  that  for  radar  systems  that  use  identical  transmitting 
and  receiving  antennas,  the  radar  cross  section  is  maximum  if 


[S]h  =  h 


(4. 36) 
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where  denotes  the  antenna  polarization  which  yields  maximum  power 
reception,  and  Xg  denotes  the  complex  eigenvalue  of  (4.36).  Moreover, 
the  scattering  matrix  can  be  diagonalized  by  a  change-of-bas i s  unitary 
transformation,  using  a  unitary  matrix  which  consists  of  the  eigen¬ 
vectors  of  (4.36)  20,42  .  The  diagonalized  form  of  [S]  is 


where  X  ,  and  X  „  are  the  eigenvalues,  with  (X  ,1  >  |X  . I .  The  maximum 
si  s2  3  '  si  1  1  s2 1 

radar  cross  section  is  given  by  [20] 


(4.37) 


The  corresponding  monostatic  power  scattering  matrix  is  thus 
diagonalized  too  : 


[P^J  =  [Sj]  [Sj]  (from  equation  (4.31)) 


By  the  invariance  properties  of  [S], 
Span  {[S]}  =  Span  ([Sj]} 

"  Ki!2+  I>s2i2 

=  Tr  {[Pd]} 


(4.38) 


Tr  { [P  ] } 


(4.  39) 
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Assuming  (4.29)  for  strict  invariance, 


Det  {[S]}  =  Oet  {[Sd]} 


=  X 


si  s2 


Therefore, 


=  Oet  { [S] } 
Span  { [  S ] } 


(4.40) 


(4.41) 


For  high  frequencies,  it  has  been  shown  in  this  thesis  that 


D  =  0.5  ej2Ar9  AF(k) 

Comparing  (4.4 1)  and  (4.42)  and  assuming  that  frequency  is 
high, 


(4.42) 

increasingly 


I* 


si 


(4.43) 


and  Arg  +  Arg  As2  2  Arg  Ap(k)  (4.44) 

Combining  (4.37),  (4.39),  (4 . 4 1 ) ,  (4.42)  and  (4.43), 

Tr  (tpn  *  ("-w) 

for  high  frequencies. 

Accordingly,  |oj  may  be  interpreted  as  the  ratio  of  the  maximum 
radar  cross  section  to  the  trace  of  the  power  scattering  matrix. 

The  power  scattering  matrix  has  been  defined  by  Graves  in  [39]  and 
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represents  the  total  power  backscattered  from  the  target  for  any 
transmitting  polarization. 

To  conclude  this  chapter,  it  is  to  be  noted  that  the  invariant 
transformation  can  of  course  be  extended  to  the  more  general  elliptic 
case.  The  scattering  ratio  has  tacitly  been  extended  in  its  definition 
to  the  general  elliptic  polarization  in  Section  4.8.  Equations  (4.7) 
and  (4.10)  which  describe  the  behavior  of  the  scattering  ratio  are 
thus  general ized. 
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CHAPTER  V 

NUMERICAL  VERIFICATION 


5.1  Data  Description 

Both  theoretical  and  experimental  data  are  available  to  verify 
the  special  case  of  the  phase-curvature  relationship  (3.12),  i.e. 
the  case  in  which  the  incident  polarization  coincides  with  one  of 
the  directions  of  the  principal  curvatures  at  specular  points  on 
the  equator  of  the  prolate  spheroid.  The  theoretical  data  was  obtained 
by  a  time-domain  synthesis  of  the  impulse  response  technique  [15]- 
fhe  solutions  generated  by  this  technique  were  checked  against 
other  theoretical  solutions  [ 8 , 1 6 ]  with  excellent  agreement  [15]. 

The  theoretical  data  were  converted  in  the  form  of  amplitudes  and 
phases  of  the  elements  of  the  scattering  matrix.  The  experimental 


data  were  measured  at  the  Electro-Science  Laboratory  of  the  Ohio 
State  University  (ESL-OSU).  The  experiments  were  conducted  [28]  on  a 


f  requency-doma i n 

SHH  (bHV  and  SVH 
spheroid  used  in 


range  yielding  the 
being  zero  in  this 
the  experiment  was 


backscattered  ret 
case).  The  'size' 
6  inch  :  12  inch. 


urns  Syv  and 
of  the  prolate 
and  the  data 


were  measured  for  two  principal  polarization  cases  in  which  a  =  0 
(VV  or  TE)  and  ot  =  r/2  (HH  or  TM)  for  aspect  angles  from  0  (nose-on) 
to  $0  degrees  (broad-side  on)  in  steps  of  15  degrees  (Figure  3.2). 
After  measurement,  the  data  were  smoothed.  Two  of  the  smoothed 


frequency-domain  data  blocks,  namely,  the  2-4  GHz  block  and  the 
4-8  GHz  block,  were  used  in  this  thesis,  covering  a  range  of  3-19 
to  12.76  in  terms  of  the  values  of  kb  (2n/\  .b).  The  error  bounds 
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on  the  experimental  data  were  specified  to  be  +  2  dtl  in  magnitude 
and  +  1U  degrees  in  phase  data. 

5.2  Direct  Verification  of  the  Phase-Curvature  Relationship 

5.2.1  Theoretical  Data 

A  typical  result  obtained  with  the  theoretical  test  data  is 
presented  in  Figure  5-1.  In  Figure  5.1(a),  the  right-hand  side  of 
(3.12)  is  plotted  against  kb  (in  steps  of  0.1  from  0.1  to  18).  From 
this  graph  It  is  clear  that  as  frequency  increases,  the  phase- 
curvature  relationship  becomes  more  accurate.  The  aspect  angle  is 
90  degrees  (broadside  incidence),  and  the  corresponding  curvature 
difference  divided  by  two  is  0.375,  to  which  the  right-hand  side 
of  (3-12)  converges.  In  Figure  5- 1  (b) ,  the  real  part  of  (1  -  ft)/( 1  +  ft) 
multiplied  with  k  approaches  zero.  Although  not  presented  in  this 
thesis,  the  real  part  itself  (without  the  factor  k)  converges  to 
zero  at  a  much  faster  rate,  particularly  at  large  values  of  k.  In 
Figure  5.1(c),  the  imaginary  part  multiplied  with  k  tends  to  the 
value  0.375.  The  scattering  chart  (imaginary  part  versus  real  part) 
is  shown  in  Figure  5- 1(d).  As  predicted,  the  plot  is  indeed  a 
spiral  which,  as  k  is  increased,  converges  to  the  point  (0,  -0.375) 
on  the  imaginary  axis. 

5.2.2  Experimental  Data 

An  extensive  amount  of  testing  of  the  phase-curvature  relation¬ 


ship  has  been  conducted  with  the  experimental  input  data.  It  was 
realized  that  because  of  the  nature  of  the  tangent  function,  a  direct 
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test  of  (3-12)  with  input  data  which  have  a  i  10  degree  error  in 
phase  measurement,  is  not  very  useful  in  a  graphical  presentation. 
Thus  the  complex  plane  plots  of  some  typical  measurement  data  are 
presented  in  Figures  5-2  to  5.**. 

The  plots  in  Figure  5.2(a)(i)  to  (iii)  are  the  experimental 
versions  of  the  theoretical  plots  in  Figures  5.1(d),  (b)  and  (c). 
Figures  5.2(a) (i)  does  indeed  hover  around  the  predicted  point  on 
the  imaginary  axis.  This  behavior  is  not  so  clearly  visible  in  Figure 
5.2(b) (i) ,  where  according  to  the  theoretical  predictions  this  plot 
should  have  given  better  results  for  the  higher  frequency  range. 

This  discrepancy  is  mainly  attributed  to  the  following  factor  : 
the  phase  error  magnified  through  the  tangent  function  gets  even 
more  magnified  through  the  multiplication  with  large  values  of  k. 

The  plots  shown  in  Figures  5.3(a)  and  (b)  are  for  the  nose-on 
incidence  case  (<t>  =  0)  for  which  there  is  no  polarization 
dependence  and  («u  -  K^)/2  vanishes.  Figure  5 -  3(b)  is  of  the  4-8 
GHz  block.  The  case  presented  in  Figures  5.4(a)  and  (b)  are  for 
$  =  45  degrees,  which  is  representative  of  a  typical  aspect.  Once 
again,  the  results  for  the  4-8  GHz  block  are  not  as  good  as  the 
2-4  GHz  block. 

It  is  to  be  noted  that  the  relative  phase  error  mentioned  in 
Chapter  IV  becomes  significant  at  high  frequencies.  Figure  5-5  shows 
the  effects  of  an  error  of  0.5  millimeters  between  the  TE  and  the 
TM  antenna  positions  along  the  direction  of  propagation.  The  error 
is  simulated  within  the  theoretical  data.  The  scattering  chart  of 
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Figure  5.1(d)  now  changes  from  a  spiral  to  a  helix  elongating  mainly 
along  the  imaginary  axis  and  away  from  the  exact  value  of  (K  -  Kv)/2. 

If  the  error  is  negative,  the  helix  elongates  in  the  opposite  direction. 
An  observation  of  the  scattering  charts  of  the  4-8  GHz  block  reveals 
that  the  plots  somewhat  look  like  a  helix  rather  than  a  spiral.  Thus, 
the  relative  phase  error  partially  explains  that,  for  the  4-8  GHz 
data,  the  scattering  chart  deviates  more  from  theoretical  prediction 
than  for  the  2-4  GHz  data.  It  i s  to  be  noted  that  in  the  4-8  GHz  block, 
there  were  many  data  points  for  which  accurate  results  were  found.  In 
general,  all  the  experimental  tests  and  the  theoretical  data  supported 
the  approximate  phase-curvature  relationship  well. 


5.3  Verification  of  the  Scattering  Ratio  D 
5.3.1  Theoretical  Data 

A  complex  plot  of  D  (imaginary  part  versus  real  part)  with  the 
theoretical  data  as  input  is  depicted  in  Figure  5.6(a).  For  high 
frequencies,  the  phase  of  0,  in  theory,  converges  to  a  constant 
value  of  2  Arg  Ap(k).  However,  a  circle  rather  than  a  cluster  of 
points  of  phases  around  2  Arg  Ap(k)  results.  This  is  due  to  the  fact 
that  the  phase  fluctuates  rapidly  at  frequencies  not  high  enough. 

Even  the  range  of  frequencies  covered  by  the  set  of  theoretical  data 
(kb  up  to  18)  is  not  sufficiently  high  to  show  the  convergence.  As 
predicted,  the  radius  of  the  circle  is  indeed  0.5-  In  Figure  5.6(b), 
only  the  ampl i tude  of  D  is  plotted  versus  kb.  Clearly,  the  amplitude 
converges  to  0-5  even  at  low  values  of  k,  in  contrast  to  the  behavior 
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of  the  phase  of  0. 

5-3-2  Experimental  Data 

Figures  5-7  to  5-12  are  the  experimental  versions  of  Figure 
5-6.  Figures  5-7  to  5-9  cover  the  frequency  range  from  2  to  4  GHz, 
with  aspect  angles  being  90,  0  and  45  degrees  respectively,  whereas 
the  rest  cover  the  frequency  range  from  4  to  8  GHz.  All  the  complex 
plots  of  D  depict  an  envelope  circle  of  radius  0-5-  The  nearer  the 
data  points  to  this  envelope,  the  more  accurate  they  are  (for  high 
va 1 ues  of  kb) . 

It  is  interesting  to  observe  that  both  Figure  5-12  and  5.4(b) 
show  that  the  data  block  of  4-8  GHz  of  45  degree  aspect  is  least 
accurate  among  the  experimental  data  blocks. 
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CHAPTER  VI 
CONCLUSION 


In  the  high  frequency  limits,  the  phase  terms  of  the  polarimetric 
scattering  data  have  been  used  to  determine  the  difference  between 
the  principal  curvatures  at  the  specular  point.  It  is  also  known 
that  the  geometric  optics  backscattered  cross-section  is  related  to 
the  geometric  mean  of  the  principal  curvatures  at  the  specular  point 
(for  smooth,  convex  objects  at  least  [6]).  Therefore,  it  should  be 
possible  to  combine  the  above  two  concepts  and  determine  the  values 
of  the  individual  principal  curvature  at  the  specular  point.  Further, 
by  judicious  use  of  differential  geometry  (e.g.  Minkowski's  formulation 
or  Chr i stoffel -Hurwi tz  formulation,  etc.),  these  curvature  values 
may  yield  the  actual  profile  of  the  scatterer. 

It  must  be  pointed  out  here  that  the  phase-curvature  relationship 
concept  does  require  a  smooth,  convex,  well-behaved  surface  structure. 
If  there  are  edge  or  other  singular  types  of  source  contributions  to 
the  backscattered  signal,  curvature  recovery  with  the  method  of  this 
thesis  will  not  b^  accurate.  However,  in  the  most  recent  development 
of  radar  target  discrimination  [5]  the  relative  phase  difference  of 
the  like  polarized  elements  has  been  found  to  provide  one  of  the 
most  important  classifiers  for  discriminat ing  a  smoothly  (undulating) 
curved  surface  from  a  d i scont i nuous I y  rough  surface  with  sharp  edges. 

It  is  found  that  there  are  several  probable  reasons  why  the  2-4 
GHz  experimental  data  give  more  accurate  results  than  the  4-8  GHz  data. 


92 


They  are  the  tangent  function  in  (3.10),  the  k  factor  and  the  relative 
phase  error,  which  all  become  more  significant  at  higher  frequencies. 
Thus,  the  phase  measurement  needs  to  be  very  accurate  and  more 
accurate  methods  of  obtaining  the  complete  broadband  polarimetric 
scattering  data  should  be  investigated. 

It  seems  that  not  only  is  the  2-4  GHz  range  valid  for  the 
first  order  correction  to  physical  optics  for  a  6"x12"  prolate  spheroid, 
but  it  is  also  a  compromise  range  between  the  high  frequency  condition 
required  by  the  curvature  recovery  model  and  the  drawback  to  lower 
frequencies  required  to  prevent  critical  magnification  of  measurement 
errors.  At  too  high  a- frequency,  the  product  of  k  and  the  tangent 
function  will  lead  to  erroneous  results  from  measurement  data. 

Although  the  curvature  recovery  mode)  is  restricted  to 
conducting  scatterers,  it  may  be  possible  to  extend  the  model  to 
dielectric  scatterers,  since  the  space-time  integral  equation,  on 
which  the  model  is  based,  has  recently  been  extended  and  applied 
to  dielectric  bodies  by  Mieras  and  Bennett  [371- 

In  concluding,  both  the  amplitude  difference  and  the  phase 
difference  of  like-polarized  elements  of  the  scattering  matrix  in 
linear  polarization  tend  to  zero  at  high  frequencies,  yet  the  phase 
difference,  however  small,  does  contain  curvature  information  of 
the  scatterer.  Regardless  of  the  type  of  orthogonal  polarization 
bases,  the  phase  sum  tends  to  a  constant  value  which  is  twice  the 
argument  of  the  Fourier  transform  of  the  silhouette  area  of  the 
target.  The  phase  sum  also  tends  to  or  equals  the  argument  of  the 
scattering  ratio  defined  in  this  thesis.  The  magnitude  of  the 
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scattering  ratio,  whose  definition  is  independent  of  whether  linear, 
circular  or  general  elliptic  polarization  is  used,  approaches  0.5 
rapidly  as  frequency  is  increased.  The  magnitude  of  the  ratio  is 
interpreted  as  the  ratio  of  the  maximum  radar  cross  section  to  the 
trace  of  the  power  scattering  matrix.  The  complex  plots  of  the  scattering 
ratio  provide  a  simple  check  on  the  accuracy  of  high  frequency 
polar imetric  measurements.  Another  curvature  recovery  equation  has 
been  derived  in  circular  polarization  basis  vector  notation.  The 
curvature  recovery  model  is  proven  to  satisfy  the  image  reconstruction 
identities  of  invariant  transformat  ion.  Finally,  the  values  of  kb 
from  3-19  to  6.38  have  been  found  to  be  most  potentially  suitable 
for  curvature  recovery  of  the  6"xl2"  prolate  spheroid  (and  probably 
targets  of  similar  size  and  shape),  provided  that  polarimetric 
measurements  can  be  improved  to  a  better  accuracy,  and  that  further 
correction  to  physical  optics  approximation  can  be  obtained. 
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APPENDIX  I 

DERIVATIONS  OF  EXPRESSIONS  FOR  PRINCIPAL  CURVATURES  FOR  PROLATE  SPHEROIDS 

Let  the  semi-major  axis  of  the  prolate  spheroid  be  C,  and  the 
semi-minor  axes  be  a  and  b,  such  that  C  >  a  =  b. 

Let  7(u,v)  be  a  curvilinear  mapping  to  represent  the  surface 
of  the  prolate  spheroid.  A  set  of  parametric  equations  are 

x  =  a  sin  u  cos  v 

y  =  b  sin  u  cos  v 

z  =  C  cos  u 

In  vector  form, 

r  -  (a  sin  u  cos  v)  i  +  (b  sin  u  sin  v)  j  +  (C  cos  u)  k 

3r 

=  (a  cos  u  cos  v,  b  cos  u  sin  v,  -C  sin  u) 

3u 

3r 

—  =  (-a  sin  u  sin  v,  b  sin  u  cos  v,  0) 

oV 

By  definition, 

3_L 

3u  3u 

2  2  2  2 
=  a  cos  u  +  C  s i n  u 

.  3r  3r 
3v  '  3v 

2  .  2 
=  a  s  i  n  u 

-y  -y 

_  _  3 _r  _3r_ 

3u  3v 

=  0 
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F  =  0  implies  that  the  u,  v  parametric  curves  are  orthogonal 
to  each  other. 

Let  n  be  the  unit  vector  normal  to  the  surface. 


r  '•  r  ,  ’  r  )  r  . 

n  =  —  x  —  /  —  x  — | 

• u  ’V  3u  tv 


2  2 
(bCs in  u  cos  v,  aCsin  u  sin  v,  absin  u  cos  u) 


2  2  4  2  2  2  4  ~ 2  2  2  2  2 

•  b  C  sin  u  cos  v  +  a  C  sin  u  sin  v+a  b  sin  u  cos  u) 


2  • 

— j  =  (-a  sin  u  cos  v,  -b  sin  u  sin  v,  -C  cos  u) 
■  u 


2  * 

— j  =*  (-a  sin  u  cos  v,  -b  sin  u  sin  v,  0) 
■>v 


-aC 


2  2  2  2  2 
■u  »(C  sin  u  +  a  cos  u> 


n  . 


2*  r  .  2 
’  r  -  aCsin  u _ 

2  *  ,.2  .22  2  , 

• v  *  C  s i n  u  +  a  cos  u } 


•v  -u 


=  (-a  cos  u  sin  v,  b  cos  u  cos  v,  0) 


,2- 

J  r  „ 
n  .  — — -  =  0 

3v3u 


By  definition, 


L  =  n  . 


32r 
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aC 


...  2  2  .2.2, 
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N  =  n  . 


g2r 

3v2 


r  •  2 

-  aCsin  u 


/,  2  2  .2.2, 

/{a  cos  u  +  C  sin  uf 


M  =  F  =  0  implies  that  the  lines  of  curvatures  are  the  u,  v 
parametric  curves  chosen.  Thus, 


K 


u  E 


-  aC 


/  2  2  r 2  .  2  , 3/2 

(a  cos  u  +  C  sin  u] 


K 

v  G 


-  C 


j,  2  2  r2.2, 

a/{a  cos  u  +  C  sin  u; 


and  the  Gaussian  curvature 


K  =  K  K 
u  v 


1 


l2,x2  y2  z2  2 
a  C  {— 2j-  +  —jf  +  — jj-} 
a  a  C 


In  terms  of  x,  y  and  z, 


Ku  "  „2 


-  aC 


C4-  2  C2  2  a2  2  ^2 

[“  *  +  ~  y  +  ~2  Z  3 
a  a  C 


K  = 
v 


C2  2  C2  2  a2  2  ' ^2 

a[— 2  x  +  —  y  +  —  z  J 

a  a  C 


Without  loss  of  generality,  consider  the  x-z  plane: 


J 
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The  point  P(x,z)  with  an  aspect  angle  <p  ,  as  shown  in  Figure 
A ( 1 ) ,  can  be  expressed  in  terms  of  <(>,  by  solving 


x  =  z  tan 


and 


2  2 
*_  +  =  1 
2  p2 

a  C 


Thus , 


aC 


z  = 


/{C^tan^  +  a~} 


aC  tan  4 

X  ~  .2.  2, 


/{C  tan  4>  +  a  } 


Consequently, 


K  -  ,  -  *C 

Ku  .  _i~ 


C  tan  <j)  +  a 

2  2  2J 
C  tanl  +  a 


V3/2 


K  = 
v 


-  C 


X  r 


/rC  tan  <{>  +  a  t 

aV  2  2  2 ; 
C  tan  <J>  +  a 


For  a  2:1  prolate  spheroid,  C  =  2a.  As  an  example,  suppose 
=  90  degrees,  hence 


K  =  - 
u 


"5a 


Kv  =  ‘  7 


K  -  K  , 
u  v  J 

77 


4 
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Since  k,  the  wave  number^  has  been  normalized  in  the 
phase-curvature  relationship,  a  is  set  to  unity  here  in  conformity. 

Thus , 


K  -  K  ys  2-7C 

_u _ v  =  0.375 

2 

for  broadside  incidence. 

Other  values  of  curvature  difference  can  be  similarly 
evaluated  once  the  aspect  angle  <j>  is  specified. 
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APPENDIX  II 
COMPUTER  PROGRAMS 


(a)  for  computing  curvature  difference  and  plotting  Figure  5.1(a) 

(b)  for  plotting  the  scattering  charts  with  theoretical  and 
experimental  data 


(c)  for  plotting  the  scattering  ratio  D 
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l)*CQS(THI))  ( 

TAU  =  SQRT(T  AU) 

H  1  =48  S( H  1  ) 

GG=G 1 *C/ ( 4*SQRT ( Ml  )) 

GG=GG/( T AU**3 ) 

P=  SQR  T(  4*A*ZT*Z  T*B *B *ET*E  T*C*  C*  AX*  AK)  ( 

01  2=  (  (  (  A*ZT  )**2  )*(  8*8— C  *C  )*«(0*ET)**2)*(4*4-C*C))**2 
01 2=01 2* ( (  ( 8*ET  )**2) *( 4*4-C*C  )♦( ( C*A<  >**2  »*( A* A- 8 *9  )  )**2 
01  2=0  12-(  (  9*ET  )  ***)*((  A*A-C*C  )**2)-2*  (  (  A*ZT>**2  )*(  (  C*  AK)**2»*  (  A*A-( 
1  e*B)*  (  B*B— C  *C  ) 

0  12=  S'IRT  (012)  *P  /  (  (  A*  B*C  )  **2  ) 

Cl  2=01  2*8/2  • 

PRINT  99.  012 

99  FORMAT!  *0*.  •  EXACT  CURVATURE  DIFFERENCE  BY  2  I S  F10.6)  i 

PRINT  5 

5  FORMAT!  'U»,  8  X ,  'NORMALIZED  K»,  30X ,  'CURVATURE  DIFFERENCE  BY  2  FR 
COM  PHASE') 

CO  6  1  =  1  .  N 

K(  I  )  =K  0  ♦  C  1-1  )  *  X  S 

AKG=(PHIH(!)  -  PHI  V(  I  )  )  /  2  • 

RHS(  I) «<( I  ) *T  AN ( ARG) 

6  CONTINUE 


102 


oo  r  i*t  ,n 

7  OR  IN  T  8*  K  (  I  )  t  RHS(1  ) 

8  F2RMAT**0*.  10X,  F9.2,  30X.  FIS. 61 

PRINT  9#  PMIOEG 

9  FORMAT* »0».  »THE  ASPECT  ANGLE  HAS  BEEN  ».  F10.6.  »  DEG*  » 
PRINT  12.  THIOEG 

12  FORMAT*  •  0*  .*  ANGLE  TH  I  *'  .F10.6)  *  OEG*  ) 

PRINT  10.  KO.  KS 

10  FORMAT*  •  0»  .*  ST  ARTI  NG  NORM  K  *  '.F10.5,  »•  STEP*  •»  F10.51 

CALL  9GNPLT* 5.CURVA) 

CALL  SCALE(K.12  .0.  180.  1  ) 

CALL  SCALE*9HS.e.O  .1  80.1  I 

CALL  AX  IS(0  .0.0  .0.  ‘NORM  K  *  . -6  •  1  2  .  0.  0 . 0  .K  (  *  81)  .  K  (  I  82  >) 

CALL  AXIS*  0.0  .0.0.  *RHS*  •  3  .8  .0  .  90  •  0  •  R  HS  (  *  8 1  l.fiHS*  182  >) 

CALL  LINE(K.RHS. 180. 1. 1 .4) 

PRINT  2221  «K( 1 81 ) .«* 182 1 . RHS* 1 81  1  .RHS* 182  > 

2221  FORMAT* *0* . 4F15.S1 
CALL  ENOPLT 

222  FORMAT*  *0*  .4F14. 61 

STOP 

END 


FCSC’A  NOTERM. ID.  EBCDIC.  SOURCE.  NO'-  I  ST. NO DECK .LOAD .NON AP .NOTE  ST 
FF*CT*  NAME  *  main  •  LINECNT  *  50 

SOURCE  STATEMENTS  *  70. PROGRAM  size  * 

NC  DIAGNOSTICS  GENERATED 


5908 
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1  A  h  u  r  A  „  E  2.9  DATE  =  8213A  10/53/3* 

^  0-AUOHU.1I‘S  SCAT  T*  2 1  NO  CHART 

<-  file  name  ::::::::::: ikchart .theory .cntl 

r  )ft.rA  ra  St  H£AO  F9 J 4  V90 . THFOP Y . OA T A  AND  H90 • THFOpY .DATA 
C  N  =  <n.  CF  FREQUENCIES  OR  DATA 

c  re  sphe  ioio 

HMEN51CN  AMPVtlTO).  AMPH(liO).  PMtV(tAO).  PH  I H  (  19  0) 

w  r  a  L  K(H2).<3.<i.  OP  (  1  3 5  )  «  IPI182).  K^d^l,  <1P(!S2) 

CJ^’LilX  CJ.  Cm*_X,  CEX3  •  PNOM.  RDFN,  F  (  1  °0  ) 

0=1  IT 
<0=0. 1 
<3*0.1 

C  J  =C  MP1_  X  (  9  »  3  «  1.0) 

.>I=*.«ATAN(  1.  ) 

THtjec=oo. 

PHt  >£G=90. 

n  _*AO(  i,  l  J  (  AMP  V  (  I  )  .PHI  V  (  I  )  .  I  *  I  »  N  ) 
tlAJl  >.  1  )  (  AXPHU  |  I  )  .  1*1  .  N> 

t  pc;aaT(0fio.6) 

30  >.  I  =  t  .  N 

■<(  1  ):<CMI-1  )  •  <  3 

p  J  [rF  =  POIH(  I  1  -  J  -1  I  V  (  I  ) 

a  >=wph(  ()/a‘«p><  r ) 

*■  NO  4=1  .-AM*CEX  >{CJ*»01FF) 

3  OF  0=1.  *\>*  *Ce<3<CJ*P0IFF) 

F  (  I  )  sON'JM/PQEN 
9  3(  l  )  S‘«i  AL  (  F  (  I  )  J 
t3(  l  ) =4  I  MAG (r  (in 
Ol)  I  »=<  (  I  )  *  TP  (  I  ) 

<  IP<  I )-< (  I  3  *  I  3 ( I  ) 

2  C  C  \  T  I  NvjC 

p-MNT  0 

j  KOUTCJ'.'NCl'AU'EO  <•  .  I  )X  ,  •PEAL'  .  1  OX  .  *  I  M  AG  I  (ADV  .  ,?5X  ,  «*oe*L' 

Cl  >X,  •  K.  I  *4  A  Ci  •  1 

DO  \  1  =  1  iN 

4  »  (I  NT  5,  G(I),  (Pin.  I^ll).  KRP(I).  KJP(I) 

■3  FC  1MATI  •  0*  ,F  l  3.  !  ,2  l  X  .  F  |  A.  7  ,  •>  X  ,  F  |  A  .  7  ,  11X  .  F  1  4 . 7 , 1  G  X  ,  F  I  4  .  ?  > 

CAUL  iO.NPL  T  (  j  «  CHART) 

‘  )<  tcIC3  SCALE  FOR  AuL  I  «  AG  I  NA  R  Y'  PF  AL  CL  J  ST  ER  PLOTS  'in*  JSE9 

c  hh£u  plots  have  n:  jn i e i e d  scale 

_  Tf.  nrvUr  ON  l  F  l  E  0  SCALE  FOR  I‘*AO-FEAL  CLUSTER  RL^TS, 

3  J  JST  )ELE  TE  THE  ^CL.CXING  A  FOPTRAN  STATEMENTS. 

AND  INSERT  d  SOOFJ.riNE  SCALE.  AS  FOLLOWS: 

C  c  all  SC  ALE ( <PP . 3 . 0 . I  SO  . \  ) 
u  v-ALL  GCAlEIK  I  P.3.  7.  130,  l  ) 

/.WP(  101  )  =  -2.  1 
<IP(  ISl  1  ) 

< * p (  m3i  - o . s 

<  I  F(  l  j,>  >  =0.5 

CALL  A/I3I0.T.A.0.'  <WE  AL  •  .  -3 . 8  .0  .0 .9  ,r  RP  (  1  S  1  )  •  K  O'1  (  IS?)  ) 

CXLl-  ax!  3(4.9.  ).0  ,  •<  I  MAS'  .-.fl.O.OO.O.KIOCISI  >.KIP(  IS?)) 

CAUU  L roF( <?P  ,<i =». 1  SO .  1 .0 • A) 

CALL  t'-JPLT 

CALL  1  'ON  PL  T  (  5  ,  <E  AL°  T  ) 

CALL  SCALE  1<  ,  1  ?.  0  .  HO  •  1  ) 

CALL  SCALE <K?P  .  i. 0  .  ISO  .  1) 

CALL  AX1  3<  0. 3,  ).9,  *NGPM  K  •  .- 6  ,  l  2 . 0 . 0 . 0  ,  K  <  I  S  l  )  .  K  (  1  S  ?  )  ) 

CALL  AXI S<0.  3,3.0.  'KfiF*  ,1,8. 0.90.0. K RP 11  81  >  ,KNP(  \S2)  )  j 

CALL  LINE(K,<  T  3  « 1  )0*  1.0*4)  J 

CALL  enoplt  1 


{ 

t 
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C  ALL 
CALL 
CALL 
CALL 
C  ALL 
CALL 
CALL 
PR  I  Nt 


HLNPL T  { 4  *  I  MAG) 

SCALE (<  .  12.0.130 . 1  ) 

SCALE(KIP*3. 0.130. 1) 

AX  I S( 0.0. 3.0 . • NORM  K*,-6.t?.0.0.n,K(tai),K(I32)  > 
AX!S( 0. 3. 1.0.  •KIM*  .3.3.0. 00 .0.KI°(1S1  )»KIP(132)  1 
LINE<  K .< I  3. ISO. 1 . 0.4  1 
END  PL  T 
553.  TMtJES 


o3  ^nq^ATCO'.  *  A  s  JL*£  ThI  =  *,F9.3,  •  DE 

PR  1  N  T  555.  PHlJfG 

>f  5  FCR AAT ( • 0* . • THE  Ai’ECT  ANGLE  =  ',F6. 

PRINT  556.  XI. <i 

*«TF  F0R4ATI  *0*  .*  STARTING  NORMALIZED  K  =• 
PRINT  991 

i‘M  FORMATI'O*.  •  ME  1RIGINAL  VERTICAL  P 
PH  I  NT  99?.  (A4RV(  11 .RHIV(  I ) . I=1,N) 

">2  E0<V4AT<  •  0*  .  IOFIO.3) 

PRINT  Q9  3 

<3  FCR4AT(*0*.  'THE  ORIGINAL  HORIZONTAL 
»01NT  0  94  .  (  »«H(  I  )  ,PHIH(  I  )  ,  1  =  1  ,N) 
FC'HATI  'T,  13F1D.1) 

GTO° 

E  NO 


DEGREES*  1 


e.  ■«.  *  DEGREES  •  ) 

-  *  «E 10.5.  *  ST;3s * . E10.5) 
POLARIZATION  CATA:*) 


°OL  AR I ZAT  JON 


c'rH'CT*  NCTLRM,  I  D.  R  JCOI  0. 3  3  JRCE  i‘i%  I  5T  ,  NODEC  K  ,  L  0  AO  .  noma  p  ,  sj  OT  F  5  T 
:r(TCT»  NAV;  ■=  MAIM  •  L  I  NEC  N  T  =  5" 

SOURCE  STATEMENTS  =  69. PROGRAM  SIZE  =  10796 

N"  DIAGNOSTICS  GEMERATiJ 


r.nnriiin  r.r>  nr. 


105 


RELEASE  2.0  MAIN  DATE  =  B2 1 3*  11/02/31 

ChAOOHUR  1  *  S  5CA1T£?ING  chart 
file  name :::::::::::: kchart.cntl 

N  =  NC.  CF  F  PE  QUE  NC  l  E  S  JR  DATA 
■’.(CUTE  SPHEROin 

DIMENSION  AMJV(’03>.  A«*>H(200),  PHIV(200>.  PHTH(ZOO) 

REAL  <(202)#  R3(232),  I°(202).  <RP(202).  MP<20?) 

COMPLEX  CJ.  CM3. <#  CFXP.  RNUM,  RDEN,  M200) 

Ns  2  >0 
C  I  N  =  3  . 

C=C I N*2 .54/1 C  3. 

3=C 

A  =  3  #*C 
A IN  =  2  #*C  IN 
F  0  =  2  . 

S=0. 01 

cj=:*«plx(o.o.  l.oi 

p 1  =  » . *A  T AN{  1  .  ) 

r hi oec=90. 

=*Hl  DEG=90. 

1  fc  A  L»  (  3  .  1  )  (  A  M  a  V  {  l  )  .  '•H  l  V  (  I  )  «  I  =  I  «  N  ) 

PEA?( 9. 1  ) ( AM3H(  I ) . OHIH<  1  > .  1  =  1  . N ) 

1  F  CRM aT {  1  OF  10.3) 

TC  2  1=1 #N 

<(  1)=<2.*PI/3.'*)*(FJA-(J-1  >  *3  )  *9 
°0  I  FF  :PH  I  r-t  (  I  )  —  3H  I  V  (  I  ) 

A«rlC.**((A'nH(l  )-»(3V(  !)  )/lO.) 

3nj4=1 . -AM*C£X3(  cj*°0  Icc ) 

'OEMal .♦AMACEX3( CJ* ’OIFF) 

r  (  1  )  sRNUM/POEN 

K°(  I  Isfl^AKF  (I)  | 

Ia(  I  ) aA I MAC( F  (  I)  ) 

«°(  I  )=<(  I  )  *  1  3  (  I  ) 

<lP(I)=K(I)*Ia(I) 

2  CONTINUE 
PRINT  3 

3  F-ir  4AT  (  •  0*  •  •  NC'MAH  ZED  K  *  .  1  0  X  .  *  R  EAL  •  .  I  RX  ,  •  I  V  A3  I  N  A  R  V  •  ,  25X  •  'K’CAL'  . 

Cl  9X.  • K I  MAG  *  ) 

no  a  i=i •  n 

A  PRINT  9.  MI).  IP( I),  IP(I),  KRPCI),  <ro(I) 

5  Fn-)4AT('0',F13.  j,21X,Fl4.7,5X,F14,7,13X.FlA.7,l5X,FlA,7) 

CALL  3GNPL  T  (  5  .  CHART) 

UNIFIED  SCALE  for  AuL  I  MAG l NARV-ofal  cluster  plots  NOX  used 
j T HER  PLOTS  HAVE  \0  UNIFIED  SCALE 

TC  REMOVE  UNIFIED  SCALE  FOR  TMAG-REAL  CLUSTER  PLOTS. 

JJST  JELL  TC  THE  FOLLOWING  4  FORTRAN  STATEMENTS, 

AND  INSERT  2  SUJRCl  TINE  SCALF,  AS  FOLLOWS: 

CALL  SCALE(<RP.R.  3.  200.1  ) 

call  scalE(<  i» »e .  o. 200 .  i ) 

<RP( 201  )=-a.O 
<i°( 201 ) =-a. o 

<7P( 202 ) =2.0 
KIP( 202 )  =  2  .0 

CALL  AXI  5( 0. 0  ,  0. ♦KPE AL*  .-5  .n  .0 .0.0 ,KR  3(201  ) »<RP( 20  ?)  ) 

CALL  AXl S ( 4  •  0  •  1.0.  • K I  MAG*  ,3.0.O,9O.O,KIP(2O1),KIP(2O2)) 

CALL  LINE(<RP.<I P.200. 1 .0. 4) 

CALL  ENOPLT 
CALL  3GNPLT  (  6  .  (cAL°TI 
CALL  SCALEIK. 12.0*200. 1 ) 

CALL  SCALE<*RP,3. 0.200, 1  ) 
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CALL  AX  I  St  0 .3 ,  3 .3, ‘NORM  < * , - 6 ,  1 2 . 0 . 0 . 0 , K < 2 0 1  ) . K C 2 0 21  ) 

CALL  AXIS! 0.0. 3.3. *KRE*  , 3  « fl. 0 , 90 . 0 , K RP l 2 3 1  >  .KRPI 20  2)  ) 

C  ALL  L INEIK ,<RP, 200. I .- 1 . 4 ) 

CALL  CNOPLT 

CALL  <3GN0L  T  I  4  «  l  MAG  ) 

CALL  SCALE  I<»  1  2.0.200.1  ) 

CALL  SCALECK  IP, 3.0, 200  .  1) 

CALL  AXI  SI  3. ‘3,  ).  3.  •  NORM  * •  ,-b , 12 .O.O.O.K <201  > .K ( ?02>  ) 

CALL  AXIS10.3.3.'D.'<IM».3.S.0.90.0.KIP<201).<IPI20?)) 

CALL  LINE<k .< IP.200. 1 .- 1 .4) 

CALL  ENOaLT 
PRINT  55  4#  CIN. AIN 

354  FOR  1AT{  •  0*  «  •  P-JCL  ATE  SPhETOHi  3*Cs*.F7.J,»  A=*.F7.3.»  INS») 

PF  I  NT  553.  TMI  )EG 

553  F  PR  N  A  T  (  •  0  *  .  • ANGLE  THI  s*.FS.3.  •  OEGPEES*) 

PRINT  555.  PHHcG 

C  6  5  FORMAT!  •  0*  .  •  The  A  i 3  fj  CT  ANGLP  =  '.rp.i.i  DFCREPS') 

PRINT  55 1> .  F0.5 

566  F  OR  naT(  *0*  »•  STARTING  FRET  = •  . F 1 0 . 5 , •  S  T£  »s •  . F  t  9 • 5 •  •  CHZ •  ) 
PRINT  09  1 

'••Cl  F0TNAT(*0'.  •  Th£  ORIGINAL  VCPT  ICAL  PCLARt^ATITN  CATAtM 

PRINT  OR2. I A*Pv(  I)  .PHI  VI  I  )  .!  =1  *  N 1 
992  FCRNATI'O*.  13-13.3) 

PRINT  993 

>93  FORAATI'O*.  *Th£  ORIGINAL  HOP  I  CENTAL  POLARIZATION  lAT*!'  ) 
PRINT  9  94  .  (  AV  J(-<  I  )  .PHI  H<  I  >  .  1=  1,  N  ) 

>94  FORMAT! *0* .  13=10.3) 

STOP 

F\n 


EFFECT*  NLTERM.  I  0.E0C3IC.  33'JRCE  «NOL  1  ST  .NCOFC  K  *  L  3  AD  •  N  n  v  a  o  ,  nOTTST 
EFFcCT*  NAME  =  MAIN  ,  LINECNT  =  50 

*  SOURCE  STATEMENTS  =  76. PROGRAM  SIZE  =  11R1C 

•  NJ  .DIAGNOSTICS  GFMF  R  A  T  e  O 


\nnc.r, 


lease  a.o 


MAIN 


DA  TF  *  3213* 


1 |/0?/0t 


ICC  *5  SCATTERING  RATIO 

name ::::::: i : : : : j rat  in .the  try .cntl 

JATA  ro  ri£  REAO  FROM  V9<3 . THEORY . DATA  ANO  H90 . THEORY. O A TA 
N  -  NO.  CF  FREQUENCIES  OR  DATA 
•  *K  "'L  A  TF  SPHEROID 

OtMENSION  AMPYI13JJ.  AMPM(|33>.  PHIV(ISO).  PHIH(!3<3) 
t'cAL  K  (  1  82  )  .  KO  «  <  S ,  RP<132>.  IP(192).  AMP<13?>.  PHASE  I  l  32  ) 

(•DUPLEX  CJ.  CMPt-X.  CEXP.  R  NU  **  •  ROEN,  R  AT  IP  (  I  30  ) 

Na|  .10 


KO* 3.1 
K?s0.l 


C  J=C  MPLKO.O.  1.0) 

P  1  =  4  .  #AT  AN  {  1  .  ) 

THI 3EG= »0. 

PHI 3EG=90. 

wEAiHa.l  MAMPVdl  .PH(V(  I)  .  f:l.N) 

REAQ1  9,  l  )(  AMPHl  I  )  ,  ->HIH<  I  )  .  1  =  1  ,  N) 
l  F ORMAT (6F1  0.6 ) 

PRINT  388 

jv*e  FTRMATI  •  0*  .  •  NC  i*  AL  I  CFO  <•  .  3X  .  •  REAL*  .  1  5  X .  •  I  M  AG  I  N  APV  •  ,  l  0  X  •  »  » HA  IF  * 
C.I4X#  *AM  PL  I TUJ= •  ) 

TO  •?  1-1 *N 

<(  I  l=KOM  I-l  )  *<  5 

r  MJ  Vs AMPH (  l  )  *CI<°(CJ*RH  |H(  I  )  )  *  AM PV  I  I  )*CE  XP I  CJ*PHIV<  I  )  ) 

<-  OEN-AMPH(  !)*«»♦•  3  M  PV  (  I  )  •*  2 


(•ahci  i )  — r mjm / t 3EN 
M’(I  )=FE  AL  I  RAT  I  J(  I  )  ) 

IP<  I  >=AI VAGI  PATIO!  I  )  ) 

A  VP  3  Qst- P  (  I  )**2  MP(  t  )**2 
A«?t  I  ) = S  OR  T (  AV3j  J) 

T  A  N= I P(  I  )/PP(  I  ) 

PmASEI  I  )  =\TAM  TAN) 

print  44Q.M  I  )  .RP(I)  ,  !P<  r  )  ,  1HAS£  (  I  ).  «“°(  I  ) 

2  C C NT  I  NOE 

V  "  y  F  c  R  A  A  T  (  •  0  •  .  F  l  0 . 3  .  \  (  5X  ,  F  l  4 . 7  )  ) 

CAUL  l)0\lPLT(5.  RATIO 
PP< 1 81 >--l .0 

iP(ian*-i.o 

T=>(  182)=<3.?5 
I  P(  1  12  >  =  0.25 

CALL  AXI  S(  0  .0  .  A.O.  •  -IE  AL  •  .-4  .S.0.0.0,  RP(  l  81  )  iPP  {  1  3?  )  ) 

C  ALL  AXlS(A,3.J.3.MMAG».*,a,0.  >0.0.  I'M  t  81  )  .  I  »  (I  3  ?  )  ) 

CALL  L  INF ( RP,  1°,  l 30.  I  .  I  .4 ) 

CALL  E  NO  PL  T 

CALL  dO.MPLTI  T  «  TAT) 

A«P( 181 ) =0.0 
A  VP(  l 82 ) =0.5 
XI 1 81 >=0.0 
K(  1 32  >=2  .5 

CALL  AXIS!  0.3.  3.0*  *  NORM  K», -8.3.0. <3.0, <(131  ).<(13?)> 

C  ALL  AX  I  5(0 .0 , 3.0,  • AMPL  I TUOE • . 9 .2 . 0 . 30 .0 , A M P ( l 8 1 ) . AM »( 1 8^ > ) 
CALL  L I NE (<  .  AMp, 130. 1*0.4) 

CALL  E  NO  PL  T 
3PINT  553,  TUJEG 

jj3  FPRMAT(*0».  *  A  a  5lE  TH!  =  ',F3.3.  •  DEGREES*  ) 

PRINT  585.  PHI  3EG 

iC.  5  FORMAT  (  •  0  *  ,  '  THE  ASPECT  ANGLE  =  *,FS.?,  '  PE  ORFF  S’) 

PRINT  556.  <  0 . <  £ 

5EC  FORMAT!  •  0 S TART ING  NOOA’ALIZEO  K  =»,F10.5,*  S  TEo=  •  .  FI  0  •  5  ) 


E 


8 


PRINT  991 

991  F  JRMAT  !  •  0*  »  •  T  H£  JRIGINAL  VERTICAL  POLARIZATION  DATA:*) 

dr i nt  992. ( am p v ( I > .PHtV! I ) . 1*1 .N> 

9vi2  FORMAT!  '  0*  •  1  OF  10.3) 

PRINT  99  J 

ic.j  F  ia«AT(  •  0*  .  *  TH  :  ORIGINAL  HORIZONTAL  POL  API  TATION  OAT  A  J  •  ) 
PRINT  994.1  AMP -l!  I)  »  PH  I  H  (  I  )  »  l  *  l  »N) 

>  >4  FORMAT!  *0*.  1  Jr  1  O •  3) 

STOP 
F  NO 
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APPENDIX  IV 
EXPERIMENTAL  DATA 


(obtained  from  Dr.  Jonathan  0.  Young, 
Electro-Sc ience  Laboratory, 

Department  of  Electrical  Engineering, 
Ohio  State  University, 

Columbus,  Ohio,  1 982 ) 
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